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SECTION  I 


INTRODUCTION 


1.1  COMPOSITE  MATERIALS 

The  idea  of  a  composite  material  is  a  very  broad  one  indeed. 

Generally  speaking,  a  composite  material  is  any  material  having  more  than 
one  distinct  phase  which  exist  together  on  a  macroscopic  scale.  This 
broad  definition  covers  materials  as  diverse  as  concrete,  fiberglass, 
and  even  the  skeletal  structures  of  animals .  In  these  examples,  two 
or  more  components  exist  together  to  create  a  "higher  performance" 
structure  than  could  be  made  of  one  of  the  components  alone.  In  this 
context,  higher  performance  signifies  a  higher  strength- to-weight  or 
strength- to-volume  ratio  of  the  material. 

An  important  class  of  composite  materials  today  are  those  in  which 
high  strength  fibers  are  molded  into  a  general  matrix  material.  This 
category  is  still  quite  broad  in  that  it  includes  such  varied  materials 
as  steel  reinforced  concrete  and  fiberglass.  This  report  is  concerned 
with  that  general  type  of  engineering  material  in  which  continuous 
nonwoven  fibers  are  oriented  in  a  matrix  in  a  very  specific  arrangement 
so  as  to  produce  a  high  structural  efficiency.  The  fibers  used  are 
typically  extremely  strong  and  possess  a  high  modulus  of  elasticity. 

Some  common  fiber  materials  axe  caxbon  and  boron,  for  example.  The 
matrix  material  ranges  from  organic  materials  to  those  of  a  ceramic  or 
metallic  nature  depending  on  the  purpose  of  the  overall  composite. 

Typically,  the  fibers  are  the  load  carrying  component  of  the 
structure,  while  the  matrix  serves  to  hold  the  fibers  in  place  and 
provide  a  means  of  transmitting  the  load  to  other  parts  of  the  material. 

As  one  example,  graphite,  which  is  among  the  most  common  fiber 
materials  in  use  today,  has  an  ultimate  strength  on  the  order  of 
500,000  psi  [-L]?  When  it  is  then  noted  that  its  density  is  one-fifth 
that  of  steel,  it  is  evident  that  a  material  made  from  it  can  have 
extremely  good  strength- to-weight  characteristics.  In  addition,  graphite 
has  a  modulus  of  elasticity  two  to  three  times  that  of  steel,  and  so  these 
composite  materials  are  extremely  stiff. 

It  is  because  of  the  good  strength- to-weight  and  stiffness  charac¬ 
teristics  that  same  of  the  first  applications  of  fibrous  composite 
materials  and  the  bulk  of  the  research  took  place,  and  still  do,  in  the 
aerospace  industry.  Aircraft  radomes  were  made  of  laminated  fiberglass 
because  of  the  high  strength  needed  and  the  desirable  characteristics 
with  respect  to  electromagnetic  wave  propagation.  In  some  applications , 
solid  rocket  motor  casings  were  made  out  of  glass  fibers.  Similar 
casings  are  now  being  designed  with  much  stronger  carbon  fibers.  In 
short,  composite  materials  are  being  used  extensively  today  throughout 


*  Number  in  brackets  refers  to  reference  number. 


the  aerospace  industries  and  are  being  chosen  more  and  more  for  the 
primary  structural  components . 

These  materials  can  technically  be  described  as  nonhomogeneous 
because  of  the  distinct  phases  present  and  nonisotropic  because  the 
overall  properties  are  dependent  on  the  direction  being  considered. 
However,  on  the  macroscopic  scale  where  the  dimensions  of  the  material 
are  large  canpared  to  those  of  the  fibers  and  the  spacings,  they  can 
be  considered  as  homogeneous.  Here,  this  means  that  properties  are 
essentially  uniform  on  a  macroscopic  scale  and  thus  do  not  vary  with 
position,  although  they  may  still  vary  with  direction. 

The  properties  within  each  component — fiber  or  matrix — may  vary 
with  direction  although  typically  they  do  not,  except  for  a  few  of  the 
fiber  mechanical  properties.  The  thermal  properties  of  the  two  com¬ 
ponents  typically  do  not  vary  with  direction,  and  none  of  the  properties 
varies  with  position  within  the  component.  Thus,  taken  seperately,  they 
can  be  considered  as  homogeneous  and  isotropic. 

Typical  composite  materials  of  the  type  being  considered  for  current 
use  are  constructed  of  many  layers,  or  laminae.  Each  lamina  is  formed 
individually  with  the  fibers  being  laid  down  in  a  parallel  fashion 
and  then  being  covered  by  the  matrix  material.  These  lamina  which  are 
one  fiber  thick  are  then  stacked  together  to  form  a  laminate,  or 
composite.  The  stacking  is  done  so  as  to  have  all  the  fibers  parallel 
and  thus  have  very  pronounced  directional  characteristics,  or  the 
layers  can  be  rotated  with  respect  to  each  other  to  give  more  isotropic 
characteristics.  The  term  quasi-isotropic  is  often  used  when  an  at¬ 
tempt  has  been  made  to  arrange  the  laminae  to  give  a  laminate  whose 
properties  are  nearly  directional-independent. 

The  simplest  composite  from  the  geometrical  standpoint  is  one  in 
which  all  the  fibers  are  parallel  and  are  stacked  in  a  square  array  as 
shown  below. 
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Figure  1.  Square  array  cross-section. 


For  the  arrangement  of  Figure  1,  it  can  be  seen  that  there  are 
three  mutually  perpendicular  planes  of  geometrical  symmetry  in  the 
material.  In  the  figure,  two  planes  are  shown  on  edge  with  the  third 
laying  parallel  to  the  paper .  These  planes  are  referred  to  as  the 
principal  planes  with  a  set  of  coordinate  axes  overlayed  on  them 
referred  to  as  the  principal,  axes. 

A  material  which  has  these  three  planes  of  synmetry  is  said  to  be 
orthotropic  on  the  macroscopic  level--a  subset  of  all  nonisotropic 
materials.  It  should  be  noted  that  a  material  can  be  homogeneous,  or 
of  uniform  structure,  and  still  have  these  directionally  dependent 
properties . 

In  the  simple  composite  structure  shown,  the  properties  such  as 
strength,  elasticity,  and  conductivity  will  be  different  along  the  three 
mutually  perpendicular  directions  since  the  two  components  will  typically 
have  quite  different  properties.  For  example,  for  the  strong  fiber  and 
weak  matrix  combination,  the  strength  will  be  much  higher  along  the 
direction  parallel  to  the  fibers  than  the  other  two  directions.  When 
the  fibers  are  isotropic,  the  overall  composite  will  be  orthotropic, 
having  the  principal  axes  with  respect  to  the  properties  the  same  as 
the  principal  axes  with  respect  to  geometry .  The  vast  majority  of 
composites  in  use  today  are  orthotropic  in  this  manner.  The  major 
exception  is  when  the  fibers  are  randomly  placed  such  that  there  are 
no  planes  of  symmetry,  and  consequently  it  is  not  orthotropic. 

A  great  deal  of  effort  has  been  extended  toward  studying  the 
mechanical  properties  of  composite  materials.  This  is  due  to  the  fact 
that  it  is  the  mechanical  properties  which  are  so  outstanding  and  which 
brought  composites  into  widespread  use.  Studies  of  thermal  properties, 
and  specifically  studies  of  heat  conduction,  are  lagging  in  the  literature 
However,  because  of  the  critical  applications  for  which  composites  are 
being  used  there  is  a  definite  need  for  studies  of  heat  conduction. 

This  report  will  present  detailed  numerical  studies  of  heat  conduction 
in  the  steady  state  in  fibrous  composite  materials. 

1.2  HEAT  CONDUCTION  IN  COMPOSITE  MATERIALS 

In  general,  the  amount  of  heat  conducted  through  an  orthotropic 
material  per  unit  area  when  in  steady  state  is  given  by  the  following 
equations  [  3] : 


In  these  equations  x,  y,  and  z  are  any  arbitrary  orthogonal  axes, 
and  the  constant  Kij  gives  the  heat  conducted  in  the  i  direction  due 
to  a  temperature  gradient  in  the  j  direction.  For  an  isotropic  material, 
only  the  values  of  Kyy,  and  Kzz  would  be  non- zero.  Hie  nine  values 
of  K  are  not  found  in  the  literature,  however,  because  of  the  infinite 
number  of  orientations  of  the  x-y-z  coordinate  system  that  can  be  used. 

Let  the  principal  axes  of  conductivity  be  denoted  by  r,  s,  and  t  so 
that  the  principal  conductivities  are  denoted  as  Krr,  Kss,  and  Kft» 

Also,  let  the  cosines  of  the  angles  between  the  principal  and  the  general 
sets  of  axes  be  denoted  by  n-y  where  i  is  the  name  of  the  general  axis 
and  j  is  the  name  of  the  principal  axis.  In  this  manner,  the  heat  flow 
per  unit  area  can  be  described  by  the  following  equations: 

■»*  *  nxrqr  +  "m’s  +  "xt’t 


«y  *  V»r  +  V»s  +  nyt«t 


1  Vr  +  "za’s  +  "zt'Jt  l2> 


The  principal  axes  have  the  unique  characteristic  that  a  temperature 
gradient  along  one  axis  will  cause  a  heat  flow  only  along  that  axis.  In 
other  words,  the  values  of  Krs,  Kr-t,  Ksr,  K^r,  and  Kts  are  zero.  This 
means  that  heat  flow  along  the  principal  axes  is  described  in  the  much 
simpler  manner  as  shown  below. 


“r  "  “Krr 

*  -k„  u 

«t  *  ‘Ktt  St  (3) 
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Because  temperature  is  a  scalar,  the  temperature  gradients  expressed 
in  these  two  coordinate  systems  are  related  by: 

32  _  „  +  n  il  +  n  *T 

dr  ~  nxr  ix  nyr  £y  nzr  dz 


32  ^  „  iT  .  *t  .  32 

da  nxa  ix  nys  dy  nzs  dz 


£T  *1  +  n  *1  +  n  IX 
dt  nxt  dx  nyt  dy  zt  dz 


Now,  by  combining  Equations  2,  3,  and  4,  the  following  relationship 
is  obtained: 

-«x  =  ^Krrnxr  +  Kssnxs  +  Kttnxt  ^  ax 

+  (Krrnxrnyr  +  Kssnxsnys  +  Kttnxtnyt>  If 
+  ^Krrnxrnzr  +  Kssnxsnzs  4  Kttnxtnzt^dz  ^ 


Equations  for  q^  and  qz  are  similar  but  are  omitted  here. 

It  is  clear  that  Equation  5  is  of  the  same  form  as  Equation  1. 
Consider  the  first  term  of  the  equation  for  qx-  Then  the  following 
relationship  holds: 


K  =  K  n*  +  K  n*  +  K.+n* 
xx  rr  xr  ss  xs  tt  xx 


(6) 


From  this  and  similar  reasoning  for  the  other  eight  conductivities, 
an  important  conclusion  can  be  drawn:  when  calculating  heat  flow  in  any 
direction  in  an  orthotropic  material,  the  only  conductivities  which  are 
needed  are  the  three  principal  conductivities. 


By  knowing  the  three  principal  conductivities  and  the  temperature 
field,  the  to tax  heat  flux  is  found  by  the  following  integration: 

Q  "  ^Area  *<n)dA  “  4r«a  *  *  ^  (7) 


Here,  q  is  the  magnitude  of  the  heat  flux  vector,  q,  and  n  is  the 
cosine  of  the  angle  between  q  and  the  normal  to  the  area,  dA-  Consider 
the  x  direction  to  be  perpendicular  to  a  plane  of  integration.  Since 
the  x-y-z  coordinate  system  is  arbitrary  anyway,  this  is  no  special 
restriction.  Then,  Equation  7  can  be  rewritten  as  follows: 


«*  ’  W  V  dA  l 

"  'Area  *xik 

*  WK«  §§  +  Kxy  if  +  g>dA  <8> 


Now  consider  the  restrictions  that  the  material  be  homogeneous 
(properties  independent  of  position)  and  that  the  temperature  gradients 
be  constant  over  the  area,  A.  Then  the  following  equation  is  obtained; 


(-K  &  - 

1  xx  ax 


K  ^  - 

xy  ay 


'Sea  55)(A) 


(9) 


Therefore,  since  these  three  conductivities  are  dependent  solely 
on  the  three  principal  conductivities  and  the  geometry,  the  total  heat 
flux  can  be  easily  found  if  the  material  is  homogeneous  and  if  the 
temperature  field  is  uniform. 

A  composite  material  can  very  often  be  considered  homogeneous 
when  the  dimension  of  temperature  resolution  is  much  larger  than  that 
of  the  fibers  or  spacings.  Then  from  that  practical  viewpoint,  the 
material  could  be  analyzed  with  proper  representation  of  an  effective 
conductivity.  If  a  uniform  temperature  gradient  were  then  imposed 
across  a  large  section  of  composite  material,  the  effective  conductivities 
could  be  used  with  Equation  9  to  determine  the  heat  flux. 

When  nonhomogeneities  with  a  conductivity  different  from  the 
matrix  material  exist  in  the  material,  the  temperature  field  will  be 
very  complex  as  it  curves  around  the  nonhomogeneities.  However  when 
the  nonhomogeneities  are  very  abundant,  the  minute  details  of  the  heat 
flux  in  the  vicinity  of  each  will  be  repetitious  and  predictable.  An 
effective  conductivity  for  the  material  would  then  take  into  account  the 
net  effect  that  the  nonhomogeneities  have  on  the  overall  heat  flux 
without  having  to  consider  the  details  of  the  temperature  field  in  the 
vicinity  of  each  inclusion. 
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Only  the  principal  effective  conductivities  are  necessary  since 
by  using  them  the  effective  conductivity  in  any  direction  can  be 
determined.  The  procedure  for  calculating  these  principal  effective 
conductivities  for  fibrous  composite  materials  will  be  discussed  in 
detail  in  the  following  sections. 

1.3  STUDIES  OF  CONDUCTION  IN  COMPOSITE  MATERIALS 

The  problem  of  heat  conduction  in  general  composite  materials  can  be 
divided  into  three  basic  groups.  The  first  of  these  groups  includes 
layered  composites  in  which  continuous  layers  exist  whose  properties  are 
different  from  adjacent  layers.  The  second  group  includes  those  materials 
where  various  shaped  particles  are  dispersed  randomly  throughout  a 
matrix  material.  The  last  group  includes  fibrous  composites  with 
specific  orientations  which  is  the  subject  of  this  report. 

Because  of  the  simplicity  of  layered  composites,  exact  methods  can 
be  employed  to  solve  for  effective  conductivities.  Their  solutions 
are  in  the  form  of  relatively  simple  expressions  and  can  be  used  as 
simplified  problems  for  more  complex  geometries.  There  has  been  a 
great  deal  of  work  concerned  with  the  effects  of  thin  films 
and  interface  resistances.  This  is  not  considered  in  this  effort. 

Materials  in  which  particles  are  distributed  randomly  through  a 
matrix  are  by  their  nature  not  suitable  to  exact  studies,  and  a 
statistical  approach  is  much  more  productive.  Hashin  and  Shtrikman  [9] 
developed  bounds  for  cases  where  only  the  volume  ratio  and  component 
conductivities  are  known  and  where  the  geometry  is  quite  arbitrary. 

These  bounds  are  quite  close  for  low  volume  ratios  but  diverge  very 
rapidly  for  higher  values.  Other  bounds  in  which  the  particles  are 
specified  to  be  spheres  or  in  which  a  well  dispersed  random  mixture  is 
specified  are  discussed  by  Hale  [ 10] .  For  dilute  suspensions  of  a 
powder  within  a  matrix,  he  has  presented  equations  to  predict  the 
thermal  conductivity.  All  of  the  work  in  this  area  points  to  the  fact 
that  with  large  volume  ratios  the  geometry  of  the  material  is  critical, 
and  accurate  predictions  of  effective  conductivity  cannot  be  made 
without  that  information. 

For  fiber-reinforced  materials,  the  Hashin  [9]  bounds  are  of  course 
applicable  since  they  allow  for  any  fiber  and  packing  geometries. 

Beran  and  Silnutzer  [11]  developed  improved  bounds  for  parallel  fiber 
composites  where  the  fiber  geometry  was  specified  although  no  attempt 
was  made  to  consider  the  packing  geometry.  Elsayed  and  McCoy  [12]  then 
developed  better  bounds  for  cases  when  the  packing  geometry  is  completely 
defined.  These  bounds  showed  that  packing  geometry  is  absolutely 
critical  for  volume  ratios  greater  than  10  percent. 

Tsou,  Chou,  and  Singh  [8]  approximated  a  fibrous  material  by 
averaging  the  properties  for  each  layer  and  then  treating  them  as  con¬ 
tinuous  homogeneous  regions  with  heat  flow  transverse  and  parallel  to 


them.  This  approximation  is  satisfactory  for  volume  ratios  of  less 
than  20  percent,  but  it  deteriorates  rapidly  for  higher  values. 


Springer  and  Tsai  [13]  refined  this  approach  by  assuming  one¬ 
dimensional  heat  flow  through  the  material.  They  were  then  able  to 
integrate  across  a  layer  of  the  material,  and  arrive  at  an  expression 
for  the  effective  thermal  conductivity.  For  a  cylindrical  fiber  and 
a  square  packing  array  with  heat  flow  transverse  to  all  fibers,  the 
equation  is  as  follows: 


*eff 


^matrix 


tan 
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(11) 


As  with  other  methods,  these  equations  are  good  for  low  volume 
ratios  but  are  inaccurate  for  high  volume  ratios. 

Behrens  [l4]  developed  a  simpler  expression  (Equation  12)  which 
incorporates  volume  ratio  and  fiber  geometry  information  but  says 
nothing  of  the  packing  geometry.  For  circular  fibers  his  equation  is 
as  follows: 


*eff  _  l  +  .... 

Kmatrix  1  -  v(f^) 

This  equation  is  also  developed  by  analogy  to  the  equations  for  ef¬ 
fective  shear  modulus  [ 1] .  A  comparison  of  these  two  equations  is  shown 
below  for  a  square  array,  (Figure  2). 

Springer  and  Tsai's  equation  gives  a  somewhat  lower  thermal  con¬ 
ductivity  than  that  of  Behrens.  They  do  however  agree  with  each  other 
within  5  percent  up  to  volume  ratios  of  50  percent. 

Because  the  effective  shear  modulus  in  a  composite  is  governed  by 
equations  analogous  to  those  for  heat  conduction,  any  results  obtained 
for  the  shear  modulus  are  directly  applicable  to  those  for  heat  con¬ 
duction.  Adams  and  Donner  [16]  used  a  finite-difference  approach  to 
determine  the  effective  shear  modulus  for  the  square  array.  Fiber- to- 
matrix  shear  modulus  ratios  and  volume  ratios  were  varied  considerably. 
The  results  they  obtained  compare  well  with  those  presented  in  the  ap¬ 
pendix  to  this  report.  They  were  however  restricted  to  the  square 
packing  array. 


Figure  2.  Comparison  of  equations  for  square  array. 

Hius,  a  number  of  methods  have  been  used  to  predict  the  effective 
conductivity  of  a  parallel  fiber  composite  material.  However,  most 
were  restricted  to  the  simplest  packing  geometry  and  made  broad 
simplifying  assumptions.  Reference  1  contains  a  good  review  of  these 
various  methods.  No  attempt  has  been  made  to  study  the  0°  -90°  case  in 
which  adjacent  layers  of  fibers  lie  at  right  angles  to  each  other. 

The  present  report  will  present  complete  parametric  studies  in  each  of 
these  areas  without  the  use  of  broad  simplifying  assumptions. 

1.4  SCOPE  OF  THE  PROJECT 

*  It  is  the  purpose  of  this  study  to  analyze  the  principal,  effective 

i(  conductivities  for  various  geometrical  arrangements  of  fiber  reinforced 

|  i  composite  materials.  The  following  basic  types  of  geometries  will  be 

'  analyzed: 

t  1.)  Unidirectional  fibers  with  all  regular  packing  patterns. 

[  j  2.)  Layers  of  fibers  perpendicular  to  adjacent  layers  with 

selected  values  of  layer  separation- to-fiber  separation 
to  be  studied. 

For  these  cases  the  following  quantities  will  be  varied  in  these 
parametric  studies: 

i.)  Ratio  of  fiber  to  matrix  conductivity, 
ii.)  Ratio  of  fiber  to  total  volume. 
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SECTION  II 


PARALLEL  FIBER  MATERIALS 


2.1  PROBLEM  DEFINITION 

The  simplest  class  of  composite  materials  is  that  in  which  a-H  the 
fibers  are  parallel  to  each  other  and  which  behaves  according  to  the 
assumptions  discussed  in  section  1.3.  The  geometry  for  these  materials 
can  be  described  in  two  dimensions,  and  when  a  uniform  temperature 
gradient  exists  across  the  fibers,  the  temperature  field  that  exists  is 
two-dimensional . 

The  configurations  of  this  class  fall  into  two  basic  categories. 
Cross-sections  of  each  are  shown  in  Figure  3*  The  first  category  is 
referred  to  as  the  rectangular  array.  The  second  is  referred  to  as 
the  staggered  array  because  every  other  layer  is  staggered  horizontally 
by  one-half  of  the  distance  between  fibers.  For  both,  the  independent 
variables  from  the  geometric  point  of  view  are  the  vertical  distance 
between  layers,  the  horizontal  distance  between  fibers  in  each  layer, 
and  the  volume  ratio  which  is  the  ratio  of  fiber  volume  to  total 
volume . 
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Rectangular  Array 

Staggered  Array 

Figure  3*  Parallel  fiber  arrays 

It  is  desired  to  find  the  heat  flux  across  a  plane  within  a 
particular  composite  material.  As  discussed  in  Section  I,  the  heat  flux 
for  an  isotropic  material  is  given  by  the  following  equation: 

Q  *  -/(K^)dA  (13) 
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Here,  n  is  the  direction  normal  to  the  plane  of  integration.  On  a  macro 
scale  a  composite  is  definitely  not  isotropic.  However,  each  component 
is  so  this  equation  is  valid  when  integrating  across  a  particular 
component . 

To  find  the  effective  conductivity  for  a  composite,  a  uniform  (one¬ 
dimensional)  temperature  gradient  on  the  macro  scale  must  be  imposed 
across  a  section  of  uniform  thickness.  The  heat  flux  across  a  plane 
in  the  material  must  then  be  determined  according  to  Equation  13 .  The 
effective  conductivity  is  then  given  by  Equation  14.  Here,  fir/An  is 
the  uniform  temperature  gradient  across  the  section  on  the  macro  scale. 


K 


eff  =  9//(AJnJ 


(14) 


It  is  therefore  necessary  to  find  the  temperature  field  in  the  material 
on  the  scale  of  the  fiber  and  spacing  size  so  as  to  perform  the  integration 
of  Equation  13. 

The  rectangular  array  will  be  considered  first  as  it  is  the  simplest 
arrangement  for  a  composite  material.  Consider  a  uniform  temperature 
gradient  (from  the  top  to  the  bottom)  to  be  imposed  over  a  section  of 
this  material.  To  do  this,  let  the  top  face  and  the  bottom  face  be  at 
constant  temperature  so  that  in  general  heat  flows  from  top  to  bottom. 

Then,  on  horizontal  lines  of  symmetry  the  temperature  will  be  constant. 

This  can  be  teen  to  be  true  because  on  lines  of  symmetry  where  the 
general  temperature  gradient  is  perpendicular  to  the  line,  effects  of 
changing  conductivity  above  exactly  cancel  the  effects  of  changing  con¬ 
ductivity  below.  Then  the  temperature  along  the  line  becomes  uniform 
as  the  overall  temperature  gradient  is  uniform. 

On  vertical  lines  of  symmetry,  lines  which  trace  the  heat  flow  through 
the  material  will  be  vertical,  in  other  words,  the  temperature  gradient 
in  the  horizontal  direction  along  vertical  lines  of  symmetry  will  be 
constant  at  zero.  This  is  the  case  because  the  effects  of  changing  con¬ 
ductivity  to  the  left  exactly  cancel  the  effects  to  the  right,  and  so 
there  is  no  net  effect  to  cause  the  flow  of  heat  to  deviate  from  a 
vertical  path.  This  of  course  assumes  that  the  material  has  large 
numbers  of  fibers.  Lines  of  symmetry  are  superimposed  on  a  section  of 
rectangular  array  below. 

Based  on  the  one-dimensional  vertical  temperature  gradient,  a 
repeating  pattern  can  be  identified  as  that  inside  the  heavy  lines. 

It  is  repeated  below  with  the  previously  discussed  temperature 
boundary  conditions  along  the  edges  indicated.  The  obvious  restrictions 
at  the  fiber-matrix  interface  are  that  the  temperature  is  the  same  in 
each  region  and  the  heat  flux  normal  to  the  interface  is  the  same  in 
each  region. 
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Figure  5-  Rectangular  array  boundary  conditions 


Therefore,  the  temperature  field  needs  to  be  found  in  this  area 
so  that  the  heat  flux  can  be  obtained  by  integrating  across  some  plane. 
The  heat  flux  will  be  the  same  as  that  for  all  the  other  repeating 
cells  in  the  material  and  so  the  effective  conductivity  found  for  this 
area  is  the  same  as  that  for  the  material  as  a  whole.  The  rest  of  this 
section  is  devoted  to  finding  the  temperature  field  given  these  boundary 
conditions  and  performing  the  subsequent  integration.  Figure  6  shows  a 
typical  set  of  constant  temperature  lines  for  a  rectangular  array  with 
the  particular  constants  given. 
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Figure  6.  Rectangular  array  constant  temperature  lines. 

The  staggered  array  is  considered  next  and  is  assumed  to  have  the 
same  one -dimensional  vertical  temperature  gradient  imposed  on  it.  The 
lines  of  symmetry  then  appear  as  shown  in  Figure  7. 


As  with  the  rectangular  array,  a  repeating  pattern  appears  as 
outlined  by  the  heavy  line.  The  temperature  field  behaves  along  these 
lines  of  symmetry  as  it  did  with  those  of  the  rectangular  array  and  so 
these  conditions  are  shown  in  the  sketch  of  the  repeating  area  in 
Figure  8. 
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Figure  7.  Staggered  array  repeating  cell. 


Figure  8.  Staggered  array  boundary  conditions. 
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As  in  the  case  of  the  rectangular  array,  the  temperature  field 
must  now  be  found  so  that  the  heat  flux  can  be  found  across  a  plane. 
Figure  9  shows  a  typical  set  of  constant  temperature  lines  for  a 
staggered  array  having  the  indicated  constants. 
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Figure  9.  Staggered  array  constant  temperature  lines. 

2.2  RECTANGULAR  ARRAY 

The  general  thermal  energy  balance  equation  for  a  solid  material 
is  the  Fourier  equation: 

pc|i=V(kVT)  (15) 

Here,  the  first  term  is  the  energy  storage  term,  and  the  second  is  the 
heat  transfer  term. 

When  the  problem  is  assumed  to  be  steady  state  with  no  heat  gen¬ 
eration  and  with  conductivity  independent  of  temperature  as  it  is  here, 
the  equation  reduces  to  that  below: 
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This  is  commonly  referred  to  as  the  Laplace  equation  and  is  valid  for 
both  the  fiber  and  matrix  regions.  When  this  is  written  out  in  polar 
coordinates,  the  following  expression  is  obtained: 


1  £*T 
r’?9 


I  *1  = 
r  *r 


(17) 


The  general  solution  to  this  two-dimensional  equation  is  that  shown 
below  where  A,  B,  C,  and  D  are  constants. 


T  =  £  rn(A„cosn6 
n=l  n 


3nsinne)  +  r”n(Cncosn0  +  Dnsinn0)  (18) 


Now,  consider  the  rectangular  array  which  is  the  simpler  of  the 
two  parallel-fiber  arrangements .  The  problem  to  be  solved  is  that 
shown  in  Figure  10.  The  temperatures  across  the  top  and  bottom  edges 
are  arbitrarily  taken  as  +1  and  0  respectively  and.  the  value  of  dT/dx 
must  equal  zero  along  the  left  and  right  sides.  Hie  horizontal  side 
is  of  length  a  and  the  vertical  of  length  b.  The  fiber  radius  is 
arbitrarily  taken  as  1. 


Figure  10.  Rectangular  array  problem 


Die  conductivities  of  the  fiber  and  matrix  are  Kf  and  res¬ 
pectively  and  the  conductivity  ratio,  Kf/Km  is  called  0.  Consider 
first  the  region  I  which  is  that  region  inside  the  fiber.  On  the  lower 
edge  where  0=0,  then  T  =  0  and  so  the  constants  Aq  and  Cq  of 
Equation  18  must  both  be  zero.  Furthur,  at  the  origin  the  temperature 
is  finite  and  so  Dn  must  be  zero  so  as  not  to  have  a  singularity  at 
that  point.  For  the  matrix  (region  II),  An  and  Cn  must  be  zero  but 
there  is  no  such  restriction  on  Dn. 

On  the  left  side  dT/dx  must  be  zero  in  both  regions.  This  is 
equivalent  to  saying  that  dT/d0  must  be  zero  at  0  =  «/2.  Since  T  con¬ 
sists  strictly  on  sine  terms,  the  derivative  will  consist  strictly  of 
cosine  terms,  and  all  coefficients  where  n  is  even  must  consequently  be 
zero  to  satisfy  the  boundary  condition.  Therefore,  the  temperature 
in  each  region  can  be  described  as  shown  below  where  the  constants 
have  been  redefined  to  avoid  confusion  between  the  two  equations. 


TI  -rJlV"*1"'"9) 

TII  -JltV"  +  =„r-n]sin(ne) 


Consider  now  the  interface  between  the  fiber  and  the  matrix.  At 
this  line  where  the  radius,  r,  was  chosen  to  be  1,  the  temperature  in 
each  region  must  be  the  same. 

at  r  =  1,  Tt  «  Ttt 


E  A  sin(ne)  =  E.,(B  +  C  )sin(ne)  (21 

n=i  n  n=i  n  n 

odd  odd 

From  this  expression,  it  can  be  seen  that 

A  =  B  +  C 
non 

for  all  n.  Therefore,  An  can.  be  eliminated  from  the  expression  for  Tj 
to  give  the  following  equations  for  temperature: 

TI  ’  rll(Bn  +  Cn>rn»in<n8)  (as 

odd 


TII  =  nIl(Bnrn  +  Cnr"n)sin(ne) 
od 


The  second  condition  that  must  be  met  at  the  interface  is  that  the 
heat  flux  normal,  to  the  interface  in  region  I  must  equal  the  heat  flux 
normal  to  the  interface  in  region  II .  This  condition  is  described  by 
the  following  equation: 


r=l 


—  V  ^TT 


r=l 


(24) 


When  Equations  22  and  23  are  differentiated  with  respect  to  r  and 
the  values  for  r  =  1  substituted  into  Equation  24,  then  Equation  26  is 
obtained. 


♦  cn>r 


n-1 


sin(ne)]  = 


odd 


.n-1 


Cnr“n"1)sin(ne)]  (25) 


Cn)sin(ne)]  (26) 


This  can  now  be  simplified  to  the  expression  of  Equation  29.  Here, 
the  ratio  is  referred  to  as 

OO  OO 

Pn£l[n(Bn  +  cn>sin(n0)3  =  nI«l><Bn  “  Cn)sin(ne)]  (27) 
oda  odd 


8(B  +  C  )  »  B 

'  n  n'  n 


(28) 


=  -B  (f— - 
n'P 


+  1 


(29) 


Using  this  equation  for  Cn,  the  equations  for  the  temperatures  can 
now  be  written  in  terms  of  one  constant  for  each  value  of  n  as  shown 
in  Equations  30  and  31* 


TI  =  ni1Bn(llB)rnsin(ne)  (30) 

odd 

TII  =  nIlBn^rn  "  (f^l)r’n>in(n9)  (31) 

odd 
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A  new  constant  En  is  now  defined  so  as  to  give  the  equations  for 
temperature  in  their  final  form. 


2Bn 

p  s  _ D__ 

en  P  +  1 


(32) 


TI  =  nf1Bnrn*in(ne)  (33) 

odd 


TII  =  -  <£-^)r-n>in(ne)  (34) 

odd 

These  are  the  general  equations  for  temperature  where  the  origin 
is  placed  at  the  center  of  the  fiber.  They  incorporate  the  boundary 
conditions  that  T  =  0  on  the  bottom  and  that  dT/dx  =  0  on  the  left  side. 

Now,  the  values  of  En  have  to  be  found.  To  do  this,  more  boundary 
conditions  are  needed.  All  possible  boundary  conditions  have  been 
utilized  in  region  I,  but  there  are  still  two  available  in  region  II. 
These  are  namely  that  the  temperature  is  constant  at  a  value  of  1  on 
the  top  edge  and  dr/dx  is  zero  on  the  right  edge.  These  conditions 
cannot  be  satisfied  exactly  at  all  points  since  polar  coordinates  are 
used  in  Equations  33  and  34.  They  can  however  be  satisfied  exactly 
at  a  finite  number  of  points  with  the  view  point  that  they  will  be 
very  closely  satisfied  at  all  other  points  along  those  edges.  Let  the 
boundary  conditions  be  satisfied  at  five,  or  M  +  1,  equally  spaced 
discrete  points  for  the  constant  temperature  condition  and  at  four,  or 
M,  equally  spaced  discrete  points  for  the  constant  dT/dx  condition  as 
shown  below.  At  the  top  right  corner  point,  both  conditions  are  to  be 
satisfied.  In  addition  dT/dx  is  already  made  to  equal  zero  in 
Equation  34  at  the  lower  right  corner  and  so  this  need  not  be  one  of 
the  discrete  points. 

The  expression  for  dT/dx  in  region  II  is  found  according  to  the 
following  procedure  using  Equation  34. 


Z  nE 
n=l  n 
odd 


C(^f^)rn  - 


+  (i^J.)r-n“l][sin(ne)co3(e)]  ♦ 
(^-f^-)r"n3Ccos(ne)=aiaLii]|  (40) 


dT 

iX 


^sinftn-l)®]  + 
JL)r"n"1sin[(n  + 


(41) 


Therefore,  Equation  4l  can  be  set  equal  to  zero  for  the  coordinates 
of  the  M  points  shown  in  Figure  11  to  give  M  equations.  Similarly, 
Equation  34  for  Tjj  can  be  set  equal  to  1  for  the  coordinates  of  the 
(M  +  1)  points  where  it  applies  to  give  another  (M  +  1)  equation.  This 
gives  a  total  of  (2M  +  1)  equations.  Now,  if  each  equation  is  carried 
out  to  include  (2M  +  1)  terms,  there  will  be  (2M  +  1)  unknowns  in  the 

(2M  +  1)  equations.  These  unknowns  are  E^,  E3,  E5, - , 

Thus,  the  equations  can  be  solved  for  these  unknowns  by  various  methods. 
They  can  then  be  substituted  back  into  the  temperature  equations  (33  and 
34)  to  give  the  temperature  at  any  point.  The  only  approximation  that 
has  been  made  is  to  treat  the  power  series  as  that  of  only  (2M  +  1) 
terms  instead  of  an  infinite  number.  The  assumption  here  is  that  the 
terms  not  included  contribute  to  a  negligible  value  compared  with  the 
sum  of  the  first  (2M  +  1)  terms.  As  was  shown  earlier,  this  temperature 
field  can  then  be  used  to  calculate  the  effective  conductivity  which 
will  be  done  in  Section  2.4. 


2.3  STAGGERED  ARRAY 

The  problem  to  be  solved  for  the  staggered  array  is  that  shown  in 
Figure  12.  For  reasons  that  will  be  discussed  below,  the  temperature 
across  the  top  will  be  taken  as  +1  and  that  across  the  bottom  as  -1. 

Consider  now  the  analysis  in  Section  2.2  for  the  rectangular  array. 
The  analysis  leading  to  the  general  equations  for  temperature  (Equations 
33  and  34)  is  valid  for  the  staggered  array  shown  except  that  one  minor 
change  must  be  made.  That  change  is  that  now  the  temperature  across  the 
bottom  edge  is  taken  as  -1  instead  of  0.  To  compensate  for  this,  the 
constant  1  must  be  subtracted  from  the  rectangular  array  equations  to 
give  the  following  equations  for  the  staggered  array: 


22 


n  ^ 


m 


i 


,T*-1 


Figure  12.  Staggered  array  problem. 


TI  =  -1 


+  e  E  rnsm(ne) 
n=l  n 
odd 


■  -1 


s,EnC'1-ra>rn  +  (H 

n=l  n  z  z 


•)r“n3sin(n8) 


Because  of  the  added  complexity  of  region  III,  a  similar  equation 
cannot  be  written  for  the  temperature  in  that  region.  However,  an 
analytical  solution  similar  to  that  used  for  the  rectangular  array  is 
still  possible  and  will  now  be  discussed. 

A  typical  plot  of  the  constant  temperature  lines  for  a  staggered 
array  with  a  value  of  p  =  4  and  a  volume  ratio  (fiber  volume/ total 
volume)  of  .442  was  shown  in  Figure  9  previously.  It  is  noted  that  a 
temperature  field  symmetry  exists  in  this  area.  The  symmetry  is  such 
that  the  temperature  field  in  the  top  half  is  a  mirror  image  of  that 
in  the  bottom  half.  Figure  13  illustrates  this  symmetry. 
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Figure  13.  staggered  array  temperature  symmetry. 


Consider  points  F  and  G  which  are  located  the  same  distance  above 
and  below  the  midplane  respectively.  They  are  also  displaced  by  equal 
amounts  horizontally  from  the  center  as  shown.  Thus,  they  are  sym¬ 
metric  to  point  0  at  the  center. 


Because  of  symmetry,  point  0  can  be  seen  to  have  a  temperature  of 
zero  because  there  are  equal  changes  in  conductivity  above  and  below 
and  to  the  left  and  right.  By  similar  reasoning,  points  F  and  G  can  be 
seen  to  have  opposite  temperatures. 


This  demonstrates  the  advantage  of  choosing  the  top  and  bottom 
temperatures  as  they  were. 


(44) 


It  can  also  be  seen  that  symmetric  points  such  as  F  and  G  have 
equal  values  of  bl/by. 
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(45) 
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This  can  be  reasoned  in  a  way  similar  to  that  used  in  looking  at  the 
temperatures.  Because  the  lines  of  constant  temperature  are  symmetric) 
the  gradient  ( bT/by)  values  must  be  the  same.  If  bT/by  varied,  then 
the  spacing  between  the  constant  temperature  lines  would  have  to  vary 
accordingly  which  would  cause  at  least  some  of  the  temperatures  to  be 
unsyrametric . 

Because  of  this  symmetry,  if  the  temperature  field  is  known  in 
one-half  of  the  area,  then  it  is  also  known  in  the  other  half.  Further, 
it  does  not  matter  which  half  of  the  area  is  studied  so  long  as  one 
of  every  pair  of  symmetric  points  is  included  in  that  half.  It  can  be 
seen  that  any  straight  line  drawn  through  the  center  point  will  separate 
the  area  into  halves  which  will  each  include  exactly  one  of  every  sym¬ 
metric  pair  of  points.  Because  points  F  and  G  are  symmetric  with  respect 
to  the  center,  it  is  impossible  to  draw  a  line  through  the  center  and 
have  both  points  in  the  same  half  of  the  area. 

Hie  solution  to  the  staggered  array  is  similar  to  that  for  the 
rectangular  array  except  that  the  discrete  points  will  now  be  located 
on  a  line  through  the  center.  At  these  points,  the  symmetry  discussed 
above  will  be  enforced  on  pairs  of  symmetric  points  along  this  line. 

This  will  provide  adequate  boundary  conditions  for  an  area  which  in¬ 
cludes  the  lower  left  fiber  to  be  solved  to  get  its  temperature  field. 

In  this  way  the  expressions  for  the  temperature  in  areas  I  and  II  are 
all  that  is  needed  to  solve  for  the  entire  temperature  field,  and  the 
complexities  of  region  III  are  now  dispensed  with. 

There  are  a  large  number  of  choices  as  to  which  line  through  the 
center  should  be  used  for  the  discrete  points.  The  only  restriction 
on  this  choice  is  that  the  line  should  not  cut  across  the  fibers 
because  no  point  may  be  used  which  exists  in  region  III. 

Consider  the  three  possible  geometries  in  which  the  volume  ratio 
is  made  as  large  as  is  physically  possible  (see  Figure  l4)  for  a  given 
ratio  of  the  sides  of  the  area.  Hie  cases  shown  are  where  the  fibers 
are  in  contact  along  vertical,  diagonal,  and  horizontal  lines  res¬ 
pectively.  The  best  lines  for  the  matching  points  are  sketched  on  the 
figures  as  dashed  lines. 

In  geometries  where  the  fibers  touch  on  the  diagonal  the  only 
matching  line  that  can  be  used  is  the  one  perpendicular  to  the  diagonal 
that  runs  from  the  lower  left  to  the  upper  right  corners.  Thus,  this 
matching  line  is  used  for  all  geometries  because  it  can  never  intersect 
the  fiber.  Then  in  cases  where  fibers  touch  along  horizontal  lines, 
additional  points  are  used  on  the  right  side.  When  fibers  touch  along 
vertical  lines,  additional  points  are  added  along  the  top.  These  points 
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Figure  14.  Staggered  array  matching  lines. 


that  are  on  the  top  or  right  side  introduce  boundary  condition  equations 
identical  to  those  that  were  derived  for  the  square  array.  For  complete 
ness,  these  equations  are  repeated  below. 

For  points  on  the  top  edge: 

TII  =  +1  *  nliEn^(^T^)rn  '  (46) 

odd 

Along  the  right  side: 

|£ll  =  0  =  n|inSn|(^-|-i)rn'1sin[(a-l)e]  + 
odd 

(^"2  1)r’"n"1sin[(n+l)6]l  (47) 


Tf  •  -1  +  +  (l^-t)r:n]sin(n..)  0*) 

odd 

Tf  '  -1  +  <“9> 
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Figure  15.  Staggered  array  matching  points. 
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(50) 


+  fi~i-a)rgn>ln(n9g>  + 

odd  1 

+  (l_|_£)rjn-)sin(n0f|  .  +2  (si) 

Everything  in  Equation  51  is  known' except -for  En*  Therefore, 
this  can  be  used  to  produce  one  equation  with  the  required  number  of 
terms  for  every  pair  of  points  used.  Further,  another  equation  is 
produced  by  using  Equation  43  to  specify  the  temperature  of  point  0  at 
the  center. 


-1 


+  nI1EnC(1_r6>rS  *  (L-f-S)ror'Jaln(n*o) 

odd 


(52) 
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Now,  the  value  of  dT/dy  must  he  made  equal  for  the  pai:;  of  points 
on  the  matching  line.  This  calls  for  an  expression  for  dT/dy  m  region 
II,  and  it  is  therefore  formulated  below. 


iT  _  £T  *£,  +  iT  *r 

dy  ae  ay  dy 


+  lv_n-l 


-n-1 


n=l  n  c 
odd 


( r  "n  ]c  os  ( n  0 ) 


=  sin(0) 


1^-  =  £cos(®) 


E  ^n  [(i-fi)rn  +  (^4-a)r-n]cos(ne)£S|lSl  + 

odJ  *• 

C(^~-^-^)rn~1  -  (1— j ^-£)r"n-1]sin(ne)sin(e))  (58) 


nEn  (^-§— ^)rn“lc°s[(n-l)03  + 
odd 


(^— 5— ^)r“n"1cosL(n+l)0]l  (59) 


Equation  59  is  the  simplest  form  for  bT/hy  in  region  II.  Now  the 
condition  of  symmetry  will  be  imposed  on  points  F  and  G  in  the  same 
manner  that  the  temperature  symmetry  was  imposed  to  give  Equation  6l 
below. 


_  *T 

ay  p  ay  g 


(60) 


n?  nEnp  2  &)r^“1cos[(n-l)9f]  +  (— ~  p)r^n~1cos[(n-H)ef]  - 

odJ 

(^-^-^)rg-1cos[(n-l)eg]  -  (^-f-^-)r“n_1cos[(n+l)03  *  0  (6i) 

Thus,  Equation  6l  can  be  used  to  produce  one  equation  with  the 
required  number  of  terms  for  every  pair  of  points  used. 

Therefore,  for  a  given  geometric  arrangement  the  matching  line  is 
drawn  as  was  discussed  earlier.  Then  discrete  points  are  placed  along 
this  line  and  either  the  top  or  right  side — whichever  is  appropriate. 

For  each  pair  of  symmetric  points,  two  equations  are  written  (51  and  6l) . 
For  the  center  point.  Equation  52  is  written.  Then  for  each  point  on 
the  top  or  right  side,  either  Equation  46  or  Equation  47  is  written. 

Each  equation  is  carried  out  to  the  same  number  of  terms  as  there  are 
equations  so  that  the  number  of  unknowns  will  equal  the  number  of  equa¬ 
tions.  The  equations  can  then  be  solved  for  the  unknowns  which  are 

E]_,  E3,  Ek,  - ,  E2m_3,  e2M-1  f°r  M  equations.  They  can  then  be  put 

back  in  tne  original  temperature  equations  v42  and  43)  to  obtain  the 
temperature  field  at  any  point.  As  with  the  rectangular  array,  the 
only  approximation  that  has  been  made  is  that  the  power  series  has  been 
reduced  to  a  finite  number  of  terms.  This  temperature  field  can  then 
be  used  to  calculate  the  effective  conductivity  in  Section  2.4. 

2.4  EFFECTIVE  CONDUCTIVITY 

Using  the  procedures  described  in  the  two  previous  Sections,  the 
temperature  field  can  be  determined  for  both  the  rectangular  and  the 
staggered  arrays .  This  temperature  field  can  then  be  used  to  obtain 
the  effective  conductivity. 

2.4.1  Rectangular  Array 


The  temperature  field  for  the  rectangular  array  is  given  by  the 
following  equations : 
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Tt  =  Z  E  rn8in(ne)  (62) 

1  n=l  ft 
odd 

TII  “  +  (L-r1)r"n3sin(ne)  (63) 

odd 

Here,  M  is  the  number  of  matching  points  that  are  used  to  satisfy  the 
boundary  conditions.  The  value  of  Tpx  is  good  for  the  entire  region 
II  for  the  rectangular  array. 

The  heat  flux  must  now  be  found  through  the  unit  area.  As  dis¬ 
cussed  previously,  it  can  be  found  by  Equation  64  by  integrating 
across  any  line  that  cuts  across  the  unit  area. 

9  =  -/(K  (64) 

Here,  n  is  the  direction  normal  to  the  line  of  integration. 

Since  the  heat  flux  will  be  the  same  across  any  line,  the  bottom 
edge  will  be  used  since  it  produces  the  simplest  expressions.  Thus, 
the  normal  direction  is  the  y  direction  and  the  analysis  proceeds  as 
follows : 


=  _/i(K  aii 
VKI  dy 


e*o 


>dr  -  'o(KII  Sy11 


e=o 


)dr  (65) 


Here  the  integration  has  been  split  up  into  the  two  regions  since 
their  conductivities  differ.  As  was  the  case  before,  "a"  is  the  length 
of  the  horizontal  side.  An  expression  for  bl/by  can  be  found  as  shown 
below. 


ai  -  ai  in  +  ii  ai 

dy  dr  dy  5?  dy 


(66) 
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(67) 


£i  e  cos. til  (68) 
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However,  the  integration  is  being  performed  along  the  line  0  =  0 
so  the  following  expressions  are  obtained: 

H  ,  1  i£|  (69) 

9  =  -Vo<r  if1!  „0,dr  -  KIi;i(r  f?I|9=0)dr  (70> 

Differentiating  the  expressions  for  temperature  in  each  region 
with  respect  to  6  gives  Equations  71  and  72. 


*  2E~1nE  r°oos(ne) 
5e  n=l  n 


*lll  -  2^-1nE  f(^4-e-)rn  ♦  (i-5-®-)r“n]cos(ne)  (72) 

46  ~  n=l  n  2 

odd 

For  0=0  the  integration  can  now  be  carried  out  to  give  Equation 
73  for  the  heat  flux. 

Q  *  (an-l)  -  (l~i)(a“n-i)] 

n=i  n  Hn=l  n  £.  i  •* 

odd  odd  (73) 

The  effective  conductivity  is  defined  by  Equation  7^  where  AT  /An 
is  the  overall  temperature  gradient  across  the  entire  unit  area  under 
study. 


q 

Area  ^ 


a(§i) 


By  nondimens ionali zing  this  expression  by  dividing  through  by  the 
matrix  conductivity,  Kn,  the  following  simplified  expression  can  be 
obtained : 


Therefore,  once  the  constants,  En,  have  been  determined  as  was 
discussed  in  Section  2.2,  they  can  immediately  be  inserted  into 
Equation  77  to  give  the  effective  conductivity. 

2.4.2  Sta^ered  Array 

Consider  now  the  staggered  array  discussed  in  Section  2.3.  In 
that  analysis  the  temperature  field  was  determined  in  only  half  of  the 
rectangular  unit  area  while  symmetry  was  used  to  give  the  field  in  the 
other  half.  The  two  basic  types  of  divisions  of  the  rectangular  area 
are  repeated  in  Figure  16. 


Figure  16.  Staggered  array  solution  regions. 

The  temperature  within  the  regions  outlined  by  heavy  lines  was 
found  directly  by  the  equations  of  2.3  which  are  given  below.  Ihe 
expression  for  Tjj  is  not  valid  outside  of  the  heavily  outlined 
regions . 


TI  *  -1  +  nlj'v'WnO 

odd 


(78) 


1  +  rJl  En^(JtaT^)rn  +  (1“r-ft)r"n]sin(ne).  (79) 


It  can  be  seen  that  as  long  as  a  is  greater  than  or  equal  to  1+5 
degrees,  the  entire  bottom  line  is  included  in  the  region  for  which 
the  equations  are  applicable.  It  can  be  seen  that  the  effective  con¬ 
ductivity  for  these  cases  can  be  analyzed  in  exactly  the  same  manner 
as  was  done  for  the  rectangular  array.  The  expressions  for  temperature 
are  the  same  except  for  a  constant,  and  this  constant  does  not  affect 
the  value  of  bT/bd  which  is  used  to  calculate  the  heat  flux  across  the 
region.  The  only  other  change  that  must  be  made  in  that  analysis  is 
that  the  overall  temperature  gradient  across  the  rectangular  unit  area, 
AT/b,  is  now  -2/b  instead  of  -l/b.  This  produces  the  following  ex¬ 
pression  for  the  effective  conductivity  for  staggered  arrays  where  a 
is  greater  than  or  equal  to  1+5  degrees: 


ir^  •  h  2?r1*n£<1-5-6>«n  - 

K«  aSSJ 


Consider  now  the  other  case  for  the  staggered  array  where  a  is 
less  than  45  degrees .  Equations  78  and  79  are  not  valid  along  the 
entire  bottom  line  and  so  they  cannot  be  integrated  along  it  to  give  a 
valid  result.  However  it  should  be  noted  that  the  center  point,  0, 
and  the  line  from  it  to  the  origin  are  always  included  in  the  region 
for  which  the  equations  are  valid.  Thus,  an  integration  can  be  performed 
along  this  line  to  five  the  total  heat  flux  across  it.  Further,  because 
of  the  symmetry  in  the  temperature  field  with  respect  to  the  center 
point,  the  heat  flux  across  the  line  from  the  center  to  the  top-right 
corner  will  be  the  same  as  that  across  the  line  segment  just  discussed. 
The  major  steps  to  this  integration  are  shown  below  where  z  represents 
the  distance  from  the  origin  to  the  center. 


Q  =  C-/*(K  f~)dr](2)  (81) 


Q  =  [Q  cos (a)  -  Q  sin (a)] (2)  (82) 

y  * 

Here,  the  integration  is  being  performed  along  a  line  of  constant 
a  while  and  refer  to  the  s  and  y  components  of  Q  respectively. 
Equation  8l  can  then  be  rewritten  as  follows: 

Q  »  [-cos(a) /*(K  j^|a)dr  -  sin«x)/*(K  f||a)dr]<2)  (83) 


These  integrations  are  performed  in  the  same  manner  as  was  done 
for  the  rectangular  array.  Thun,  the  following  expression  for  Q  is 
obtained: 
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The  effective  conductivity  can  now  be  found  in  the  same  manner-  as 
was  done  before. 
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Ey  inserting  the  previously  obtained  values  of  En  into  this 
equation,  the  effective  conductivity  can  now  be  found.  This  is  of  a 
similar  form  as  that  obtained  before,  although  now  the  angle  a  enters 
into  the  expression  where  as  before  it  was  made  equal  to  zero. 

2.5  GEOMETRIC  RESTRICTIONS 

The  solution  method  evolved  in  this  investigation  for  the  stag¬ 
gered  and  the  rectangular  arrays  was  devised  in  such  a  way  that  the 
effective  conductivity  could  be  calculated  for  any  physical  arrangement. 
The  equations  that  were  derived  will  of  course  give  numerical  answers 
no  matter  what  volume  ratio  (fiber  volume/total  volume  is  used.  It  is 
therefore  desirable  to  know  the  maximum  physical  volume  ratio  for  each 
of  these  arrays  for  various  angles  a. 

For  the  rectangular  array,  the  maximum  volume  ratio  occurs  when 
the  radius  of  the  fiber  is  made  equal  to  either  of  the  sides  of  the 
unit  area.  The  two  cases  are  shown  below. 


For  the  first  case,  where  a  is  greater  than  45  degrees,  the 
critical  volume  ratio  is  given  by  Equation  88.  For  the  case  where  a 
is  less  than  45  degrees,  Equation  89  is  the  valid  equation.  It  should 
be  noted  that  at  45  degrees  the  two  expressions  are  identical. 
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Figure  17-  Rectangular  array  maximum  volume  ratios. 
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critical  4 


a  <  45° 


(89) 


For  the  staggered  array,  three  different  situations  occur  which 
physically  limit  the  value  of  the  volume  ratio.  The  limiting  conditions 
are  shown  below. 


Figure  18.  Staggered  array  maximum  volume  ratios. 

The  first  case  occurs  for  angles  greater  than  60  degrees  when  the 
fibers  are  in  contact  along  a  horizontal  line.  If  the  volume  ratio 
were  made  any  greater,  the  fibers  would  overlap,  and  this  of  course  is 
impossible.  The  second  case  occurs  for  angles  between  30  and  60  degrees 
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and  is  characterized  by  contact  along  diagonal  lines.  Finally,  the 
third  case  shows  contact  along  vertical  lines  and  occurs  for  angles 
less  than  30  degrees.  The  expressions  for  these  critical  volume  ratios 
are  given  below. 
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critical 


)  v 

2tan(a) 


a  >  60®  (90) 
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(91) 
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a  <  30®  (92) 


Table  1  is  a  short  table  of  the  values  of  critical  volume  ratio. 

It  should  be  noted  that  there  is  a  large  variation  in  these  values 
demonstrating  that  some  materials  are  much  more  efficient  than  others 
in  terms  of  the  quantity  of  fiber  that  can  be  included  in  the  composite. 


TABLE  1. 

CRITICAL  VOLUME  RATIOS 

a 

Critical  Volume  Ratio 

degrees 

Staggered 

Rectangular 

25 

.733 

.366 

30 

•907 

•453 

35 

.836 

•550 

4o 

.798 

•659 

45 

.785 

.785 

50 

.798 

.659 

55 

.836 

.550 

60 

.907 

.453 

65 

.733 

.366 

2 .6  RESULTS 

The  results  for  the  rectangular  arrays  and  the  staggered  arrays 
are  presented  in  graphical  form  in  Appendix  A  to  this  report.  For  both 
cases  the  value  of  Kfiter/Kmatrix  was  varied  from  0.1  to  100.  The 
volume  ratio  was  varied  from  0.3  through  the  highest  possible  value. 
Finally,  the  fiber  angle,  a,  was  varied  in  5  degree  increments  from  25 
degrees  through  65  degrees. 
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The  accuracy  of  the  data  used  to  produce  these  graphs  was  good  to 
three  significant  figures.  This  was  deemed  to  be  quite  satisfactory 
from  an  engineering  point  of  view. 

Figures  A-l  through  A-7  for  the  rectangular  array  and  A- 21  through 
A-29  for  the  staggered  array  represent  complete  parametric  studies. 
Equation  12  has  been  plotted  on  these  graphs  as  this  equation  is  by 
far  the  most  common  equation  found  in  the  literature  for  predicting  ef¬ 
fective  conductivities.  This  equation  does  not  take  into  account  the 
varying  packing  geometry  and  so  it  is  independent  of  OC.  Thus,  only  one 
curve  for  this  equation  is  shown  on  each  graph. 

Equation  12  gives  accurate  effective  conductivities  for  conductivity 
ratios  below  5  or  for  volume  ratios  below  0.2.  For  values  greater  than 
these  it  is  obvious  it  cannot  be  used  because  of  the  critical  role  that 
the  angle  a  plays  in  the  effective  conductivity. 

It  is  evident  that  at  high  conductivity  ratios  changes  in  the 
packing  geometry  (the  angle  a)  can  change  the  effective  conductivity 
by  an  order  of  magnitude .  This  reflects  the  fact  that  the  upper  and 
lower  bounds  for  the  conditions  [9>  H,  12]  diverge  greatly.  More 
specifically,  it  can  be  seen  that  it  is  the  spacings  between  fibers 
which  are  critically  important. 

For  the  high  conductivity  ratios,  almost  all  the  resistance  to  the 
heat  flow  is  in  the  matrix  material.  Thus,  the  length  of  the  path 
through  which  it  must  flow  in  the  matrix  is  the  most  important  factor. 

As  an  example,  refer  to  Figure  A-2  which  is  for  the  rectangular  array 
with  a  volume  ratio  of  0.4.  Hie  30  degree  and  the  35  degree  geometries 
are  shown  below  drawn  to  scale. 


Figure  19*  Rectangular  arrays  for  volume  ratio  of  0.4. 
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For  both  cases,  most  of  the  heat  is  flowing  vertically  through  the 
fibers.  But  the  vertical  distance  between  fibers  for  the  35  degree 
case  is  three  times  the  distance  for  the  30  degree  case  and  so  there 
is  a  large  difference  in  the  effective  conductivities.  It  should  be 
realized  that  a  volume  ratio  of  0.40  for  the  30  degree  case  is  very 
near  the  maximum  possible  value  of  0.45.  This  fact  alone  indicates 
that  the  fibers  are  almost  touching  and  that  a  high  effective  con¬ 
ductivity  should  be  expected. 


For  the  rectangular  array,  the  curves  for  the  effective  con¬ 
ductivity  line  up  in  order  of  increasing  angle,  a,  whereas  this  does 
not  occur  for  the  staggered  array.  This  can  be  explained  by  considering 
which  fiber  spacings  affect  the  effective  conductivity  for  each  case. 
With  the  rectangular  array,  heat  is  flowing  almost  entirely  across  the 
vertical,  spacings  between  fibers.  Seme  heat  will  be  flowing  across 
the  diagonal  spacings  but  they  are  a  great  deal  larger  and  thus  a 
great  deal  more  resistive  so  the  heat  flow  across  these  spacings  will 
be  a  small  fraction  of  the  total.  Urns,  changes  in  these  spacings 
will  affect  the  total  heat  flow  very  little.  Then  because  the  vertical 
spacings  are  strictly  increasing  with  increasing  a,  the  thermal  con¬ 
ductivity  will  be  strictly  decreasing  as  is  evident  from  the  graphs. 

A  list  of  spacings  normalized  to  the  fiber  radius  for  an  arbitrary 
volume  ratio  of  0.3  for  both  rectangular  and  staggered  arrays  is  shown 
below.  These  are  the  distances  between  fiber  edges  and  not  those 
between  fiber  centers. 

TABLE  2.  FIBER  SPACINGS  NORMALIZED  TO  FIBER  RADIUS 


Rectangular  Array  Staggered  Array 


cP 

diag. 

hor. 

vert. 

diag. 

25 

2.74 

.21 

3.23 

4.70 

1.12 

1.70 

30 

2.26 

.46 

2.92 

4.02 

1.48 

1.48 

35 

1.87 

.70 

2.72 

3.47 

1.83 

1.34 

4o 

1.53 

.96 

2.61 

3.00 

2.19 

1.26 

45 

1.24 

1.24 

2.58 

2.58 

2.58 

1.24 

50 

•96 

1.53 

2.61 

3-00 

2.19 

1.26 

55 

.70 

1.87 

2.72 

3.47 

1.83 

1.34 

60 

.46 

2.26 

2.92 

4.02 

1.48 

1.48 

65 

.21 

2.74 

3-23 

4.70 

1.12 

1.70 

For  both 

arrays  the 

horizontal 

spacings  are 

irrelevant 

since  they 

are  perpendicular  to  the  heat  flow.  For  the  rectangular  array  the 
diagonal  spacings  are  much  larger  than  the  vertical  spacings  and  so 
they  are  relatively  insignificant.  However,  for  the  staggered  array 
the  vertical  diagonal  spacings  are  of  the  same  order  of  magnitude  and 
so  both  influence  the  effective  conductivity.  This  is  the  reason  why 


the  staggered  array  curves  do  not  stack,  up  according  to  increasing  Of. 

At  different  times  different  spacings  carry  most  of  the  heat  flow. 
Whenever  either  of  the  spacings  becomes  very  small  (fibers  almost 
touching)  the  effective  conductivity  will  become  large  no  matter  what 
the  other  spacing  is  like.  Thus  the  25-degree  curve  lies  above  the 
30-degree  curve  because  of  the  smaller  vertical  spacing,  but  the  45- 
degree  curve  lies  above  the  30  degree  curve  because  of  a  smaller 
diagonal  spacing. 

It  can  be  seen  that  in  general  the  rectangular  arrays  produce 
higher  effective  conductivities  than  do  the  staggered  arrays  for  the  same 
volume  ratios  and  conductivity  ratios.  This  is  the  case  because  the 
fibers  are  stacked  in  rows  parallel  to  the  heat  flow  which  naturally 
tends  to  give  small  fiber  spacings  in  that  direction.  Thus,  these 
materials  provide  low  resistance  paths  for  the  heat  flow.  The  staggered 
arrays  however  tend  to  have  larger  fiber  spacings  and  thus  lower  ef¬ 
fective  conductivities.  There  is  some  overlap  between  the  region 
covered  by  the  staggered  array  and  that  covered  by  the  equivalent  rec¬ 
tangular  array. 

For  a  given  type  of  array  and  a  given  volume  ratio,  consider  the 
various  pairs  of  angles  which  are  complementary  to  each  other  such  as 
35  and  55,  25  and  65,  and  so  forth.  For  a  staggered  array  with  a 
volume  ratio  of  0.2,  the  25-degree  and  65-degree  cases  are  shown  below 
in  Figure  20. 


o  o 


25°  65® 

Figure  20.  Complementary  geometries  for  staggered  array. 
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It  is  evident  that  these  are  really  the  same  physical  geometry 
which  has  merely  been  rotated  by  90  degrees.  Ulus,  the  two  effective 
conductivities  which  are  given  in  the  graphs  for  these  values  are  in 
fact  two  of  the  three  principal  conductivities  for  the  same  material. 
The  third  principal  conductivity  is  that  for  which  heat  is  flowing  into 
the  paper  parallel  to  the  fibers.  That  effective  conductivity  is 
merely  the  weighted  average  of  the  individual  conductivities  based  on 
the  volume  ratio  and  is  given  by  the  following  equation: 


K  «  K,., 

eff  ■  =  v  ^b.g£-  +  1  -  V 
K  K 

matrix  matrix 


(93) 


Figures  A-8  through  A-20  for  the  rectangular  array  and  A- 30 
through  A- 56  for  the  staggered  array  present  the  three  principal  con¬ 
ductivities  together  for  various  materials.  It  should  be  noted  that 
the  effective  conductivity  parallel  to  the  fibers  for  high  conductivity 
ratios  is  always  greater  than  both  of  the  transverse  conductivities 
because  the  high  conductivity  fibers  provide  an  excellent  heat  con¬ 
duction  path.  Further,  the  transverse  conductivity  curves  level  off 
at  high  conductivity  ratios  because  most  of  the  resistance  is  then  in 
the  matrix,  and  increasing  the  fiber  conductivity  does  little  to  de¬ 
crease  the  total  resistance.  However,  for  flow  parallel  to  the  fibers 
the  curves  become  more  and  more  linear  at  high  conductivity  ratios 
because  more  and  more  of  the  heat  is  flowing  strictly  through  the  fibers 
and  the  only  resistance  it  sees  is  the  fiber  resistance. 

It  may  be  noted  that  the  45-degree  curves  for  the  staggered  arrays 
are  identical  to  the  45-degree  curves  for  the  rectangular  arrays.  The 
two  cases  are  shown  in  Figure  21. 


These  are  really  the  same  configuration  with  the  heat  flow 
directions  differing  by  an  angle  of  45  degrees.  Further,  it  should  be 
realized  that  the  horizontal  effective  conductivities  of  these  materials 
are  the  same  due  to  geometric  symmetry.  3y  the  equations  of  Section  I, 


the  reasoning  can  be  extended  to  state  that  the  effective  conductivities 
of  these  materials  will  be  the  same  for  heat  flow  in  any  direction 
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CHAPTER  III 


0°-90°  FIBER  MATERIALS 
3-1  PROBLEM  DESCRIPTION 
3.1.1  General  Discussion 


Composite  materials  in  which  the  fibers  are  not  all  parallel  will 
now  be  considered.  In  particular,  that  group  in  which  adjacent  layers 
of  fibers  are  laid  down  at  90-degree  angles  to  each  other  will  be 
studied.  A  sketch  of  a  typical  arrangement  in  this  group  is  shown  in 
Figure  22. 

Fran  a  geometric  point  of  view,  the  independent  variables  are  the 
horizontal  spacing  between  fibers,  the  vertical  spacing  between  layers, 
and  the  volume  ratio.  As  was  the  case  with  the  parallel  fiber  arrange¬ 
ment,  only  those  arrangements  in  which  the  horizontal  spacing  of  al 1 
layers  is  the  same  will  be  studied.  Since  the  whole  geometry  can  be 
scaled  up  or  down  without  changing  the  effective  conductivity,  only  the 
ratio  of  vertical  spacing  (between  layers)  to  the  horizontal  spacing 
(between  fibers  in  a  layer)  is  needed  in  addition  to  the  volume  ratio 
to  describe  the  geometry.  This  ratio  will  henceforth  be  referred  to 
as  the  aspect  ratio . 

The  last  variable  which  is  needed  is  the  ratio  of  fiber  conductivity 
to  the  matrix  conductivity.  Thus ,  when  these  three  variables  are 
specified,  there  exists  a  unique  value  of  the  ratio  of  effective  con¬ 
ductivity  to  matrix  conductivity  which  can  be  assigned  to  the  material. 
That  effective  conductivity  ratio  will  then  be  valid  for  any  geometrically 
similar  material  which  has  the  same  ratio  of  fiber  to  matrix  con¬ 
ductivity. 

As  with  the  parallel  fiber  arrangements,  it  is  first  necessary  to 
find  the  temperature  field  in  the  material  so  that  it  can  be  integrated 
across  a  plane  to  give  the  net  heat  flow  through  the  material.  Hie  ef¬ 
fective  conductivity  can  then  be  seen  to  be  given  by  Equation  94  where 
AT /An  refers  to  the  uniform  temperature  gradient  which  exists  over  a 
large  section  of  the  material. 


Kef  f 


(94) 
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Figure  22.  0°-90°  Fiber  material. 

3.1.2  Heat  Flow  Perpendicular  to  Both  Fibers 

First  the  case  where  heat  flow  is  transverse  to  both  directions  of 
fibers  will  be  considered.  As  with  the  parallel  fiber  arrays,  a 
repeating  unit  cell  can  be  isolated  in  which  the  temperature  field  is 
identical  to,  or  a  mirror  image  of  that  of  every  other  unit  cell.  A 
typical  unit  cell  is  shown  in  Figure  23  where  heat  flow  is  along  the 
y-axis. 

Note  that  the  two  horizontal  spacings  are  equal  as  this  analysis 
is  restricted  to  that  case.  Since  the  aspect  ratio  was  defined  as 
vertical  spacing  divided  by  horizontal  spacing,  and  the  horizontal 
spacing  between  fibers  is  2(a),  the  aspect  ratio  is  the  quantity  b/(2a). 

For  a  one-dimensional  uniform  temperature  field,  the  top  and  bottom 
surfaces  will  be  at  a  constant  temperature.  This  is  due  to  the  fact 
that  the  geometric  symmetry  would  not  permit  the  temperature  to  vary 
on  these  planes  since  a  conductivity  change  above  the  plane  is  exactly 
balanced  by  a  conductivity  change  below  the  plane.  By  a  similar 
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Figure  23.  o°-90°  Array  problem 

v. transverse- transverse) . 

reasoning,  the  four  vertical  planes  of  symmetry  will  be  planes  of  zero 
normal  heat  flux.  Thus,  the  boundary  conditions  on  these  planes  are 
such  that  they  can  be  considered  perfectly  insulated  planes. 

The  obvious  conditions  which  must  be  satisfied  in  the  interior  of 
this  unit  cell  are  that  at  the  fiber-matrix  interfaces  the  temperature 
in  each  region  must  be  the  same  and  the  heat  flux  normal  to  the  inter¬ 
face  must  be  the  same.  Therefore,  a  well  defined  boundary  value  problem 
has  been  presented  such  that  a  unique  temperature  field  exists.  The 
solution  to  this  problem  is  presented  in  Section  3*2. 

3.1.3  Heat  Flow  Perpendicular  to  One  Fiber 

Consider  now  the  case  where  the  one -dimensional  uniform  temperature 
field  is  rotated  90  degrees  so  that  heat  flow  is  transverse  to  one 
layer  of  fibers  and  axial  to  the  adjacent  layers.  The  same  repeating 
unit  cell  can  be  isolated  although  the  boundary  conditions  are  now 
different.  This  cell  is  shown  in  Figure  24.  It  should  be  noted  that 
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Figure  24.  0°-90°  Array  problem 

■  axial-transverse) . 


in  this  figure  heat  flow  is  along  the  z-axis.  However,  because  the 
horizontal  s pacings  are  the  same,  the  case  with  heat  flow  in  the  x 
direction  is  a  similar  problem  and  would  therefore  have  the  same  ef¬ 
fective  conductivity. 

The  aspect  ratio  is  defined  as  before  and  is  equal  to  the  quantity 
b/,2a).  By  similar  reasoning  as  was  done  for  the  previous  case,  the 
boundary  conditions  on  this  unit  cell  can  be  described.  On  the  front 
and  back  planes  which  are  normal  to  the  heat  flow,  a  constant 
temperature  condition  exists.  On  the  other  four  sides  there  is  no 
heat  flow  normal  to  the  surface  and  they  can  be  considered  as  planes 
having  perfect  insulation  on  their  surfaces.  The  interface  conditions 
in  the  interior  are  of  course  the  came  as  for  the  other  cases. 

Thus,  this  constitutes  another  well  defined  boundary  value  problem 
which  can  now  bo  solved  for  the  unique  temperature  field  and  from 
which  the  heat  flux  can  subsequently  be  determined.  Then,  once  the 
heat  flux  is  known,  the  effective  conductivity  can  be  readily  obtained. 
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3.2  GENERAL  EQUATIONS 


The  general  thermal  energy  balance  equation  for  a  solid  material 
is  the  Fourier  equation: 

PC^=V(KVT)  (95) 

For  the  problems  being  considered.  Equation  95  simplifies  to: 

v*T  =  0  (96) 

This  is  commonly  referred  to  as  the  Laplae’  equation  and  is  valid 
for  both  the  fiber  and  the  matrix  regions.  This  equation  can  be  re¬ 
written  in  the  following  more  conventional  manner: 


3x  3y  3z 


(97) 


In  the  previous  cases  where  all  the  fibers  were  parallel  to  each 
other,  a  general  solution  was  written  in  the  form  of  a  power  series. 

The  constants  of  the  power  series  were  then  found  by  use  of  the  boundary 
conditions  and  the  temperature  field  was  aiproximated  by  a  finite 
number  of  terms. 


In  the  problem  being  considered  in  this  section,  the  boundary 
conditions  are  more  complicated  and  there  are  two  fiber-matrix  interfaces 
that  must  be  accounted  for.  This  added  complexity  of  the  problem  and 
made  it  impossible  to  satisfy  all  the  conditions  of  the  problem  by  a 
finite  number  of  discrete  points.  Thus,  an  entirely  different  method 
of  solution  must  be  employed. 

When  it  is  not  possible  to  satisfy  all  the  boundary  conditions  at 
once  by  Equation  97,  one  must  turn  to  a  numerical  method  in  which  the 
original  problem  is  broken  up  into  a  large  number  of  smaller  problems. 
Consider  here  that  the  large  unit  cell  is  broken  up  into  thousands  of 
smaller  cells.  Then  these  cells  will  have  very  simple  boundary  con¬ 
ditions  which  are  imposed  on  them,  and  Equation  97  can  be  written  for 
each  one.  If  enough  cells  are  used,  then  they  can  be  considered  to 
have  constant  temperatures  throughout,  and  the  temperatures  can  be 
solved  for.  The  boundary  conditions  for  the  cell  in  question  are  the 
temperatures  of  the  bordering  cells.  Because  of  this  arrangement,  the 
temperatures  of  all  the  cells  must  be  solved  for  simultaneously  and  so 
an  iterative  solution  is  called  for.  This  method  is  commonly  referred 
to  as  the  finite-difference  method. 


Consider  a  rectangular  grid  to  be  superimposed  on  the  unit  cell. 

It  is  layed  down  so  that  the  intersection  points  lie  at  the  center  of 
the  small  cells,  and  thus  the  temperature  at  the  center  is  the  temperature 
of  the  entire  cell.  Thus ,  the  Laplace  equation  must  be  satisfied  for 
each  of  these  points  in  finite-difference  form.  Figure  25  shows  an 
arbitrary  point  as  it  is  surrounded  by  its  six  nearest  neighbors  with 
subscripting  included  to  identify  the  points. 


i. j+l.k 


i 

i 


i. j-l,k 

Figure  25.  Interior  grid  point. 

The  first  derivative  of  the  temperature  with  respect  to  distance  can 
now  be  written  in  terms  of  the  finite  difference  in  temperature  between 
adjacent  points  as  follows.  Here,  only  the  x  direction  will  be  con¬ 
sidered  although  the  other  directions  are  of  a  similar  form.  At  the 
position  i  +  ^,J,k  the  first  derivative  can  be  approximated  as: 

.  .  .  =  Ti+t .  ,i.k  ~  Ti.,i.k 

dx  i+t, j,k  2 ax 

At  the  point  i  -  ^,J,k  the  first  derivative  is: 


(98) 


The  second  derivative  at  the  point  i,j,k  can  now  be  approximated 
using  the  above  expressions  as  follows: 


St  _  St 

cix  .  1  .  ,  ^X  .  1  .  . 


(100) 


i,J>k 


dx*  . 


=  Tl+l..i.k  4  Vt.i.k  ~  2Vi.v 


i.  j.k 


(101) 


Equation  101  is  commonly  referred  to  as  a  central-difference 
equation.  Ey  combining  the  similar  equations  for  the  y  and  z  directions, 
the  Laplace  equation  can  now  be  written  in  finite  difference  form  to  give 
the  following  expression: 

*  Ti.-l,.i,k'2Ti,.i,k  +  TLJ+l,k  *  Vl- l.k  ~  2Ti,.i.k 

2  Ax  2 ny 


T.  . 

_  l.  J. 


+  Ti 


-  2T. 


^=0  (102) 


If  Ax  ,  Ay  ,  Az  are  all  taken  to  be  equal,  then  Equation  103 
below  is  obtained. 


j.* 


P .  .+T.  +  T  .  4-T  4- 

i+1* j«k  i-l, j,k  i,j+l,k  j-l,k  + 


Ti,  j,k+l  +  Ti.  j,k-l)/6 


(103) 


Thus,  the  temperature  of  the  point  in  question  (point  i,j,k)  is 
given  in  terms  of  the  temperatures  of  the  neighboring  points.  Further, 
for  the  equal  spacing  here,  it  is  merely  the  average  of  the  temperature 
of  these  points.  This  is  the  finite-difference  representation  of  the 
solution  to  the  Laplace  equation  for  these  points. 

Using  Equation  103,  the  temperature  of  each  joint  in  the  regular 
grid  can  be  estimated  using  the  estimated  temperatures  of  neighboring 
points.  By  successively  re-estimating  the  temperatures  of  all  the 
points  and  by  doing  this  repeatedly,  the  temperature  values  assigned 
to  each  point  will  converge  to  values  such  that  Equation  103  is  as 
closely  satisfied  as  is  desired  at  every  point  simultaneously.  Depend¬ 
ing  on  the  number  of  points  chosen,  the  resulting  values  will  then 
approximate  the  actual  temperature  at  those  points  in  the  real  material. 
A  numerical  integration  can  then  be  performed  using  these  values  across 


any  plane  through  the  material  to  determine  the  heat  flux  and  thus  the 
effective  conductivity. 

For  the  points  just  described  which  are  in  a  uniform  section  of 
the  material,  the  Laplace  equation  can  be  written  for  each  one  by  using 
the  boundary  conditions  which  are  the  temperatures  of  the  neighboring 
points.  Thus,  a  boundary  value  problem  has  been  defined  with  very 
simple  boundary  conditions,  and  a  very  simple  solution  has  been  obtained. 
There  are  however  many  points  to  be  considered  which  are  affected  by  the 
boundary  conditions  that  are  impressed  on  the  large  unit  cell  itself. 
These  points  exist  at  the  sides  of  the  unit  cell  and  next  to  the  fiber- 
matrix  interfaces.  They  must  reflect  those  boundary  conditions,  and 
by  doing  so  the  effect  of  those  conditions  will  be  impressed  on  the 
problem  as  a  whole  and  it  will  be  forced  to  converge  to  the  unique 
solution  for  the  temperature  field  which  exists.  The  finite-difference 
representations  of  the  solutions  to  the  Laplace  equation  for  these 
various  types  of  boundary  points  are  presented  in  Section  3*3* 

3-3  BOUNDARY  CONDITION  EQUATIONS 

'There  are  basically  three  types  of  boundary  conditions  which  must 
be  impressed  on  the  various  points  that  are  being  considered  in  the 
finite-difference  representation  of  this  problem.  They  are  the  fixed 
temperature  condition  on  two  surfaces  of  the  unit  cell,  the  insulated 
condition  of  four  surfaces,  and  the  conditions  at  the  fiber-matrix 
interfaces. 

3.3.1  Fixed  Temperature  foundary  Points 

The  first  condition  of  constant  temperature  is  very  easily  handled 
by  the  finite-difference  method.  It  is  merely  necessary  to  assign  the 
desired  temperature  to  these  points  and  never  include  them  in  the  set 
of  points  which  will  be  used  in  the  iteration  scheme  and  made  to 
satisfy  the  Laplace  equation. 

3.3.2  Insulated  Boundary  Points 

Consider  now  an  arrangement  of  these  nodal  points  in  the  same 
regular  grid  such  that  Ax  =  Ay  =  A  z  but  located  at  an  insulated 
boundary  (see  figure  26).  The  position  i  +  l,j,k  is  where  a  point 
would  exist  were  it  not  for  the  fact  that  the  insulated  boundary  is 
the  edge  of  the  unit  cell  and  no  nodal  points  exist  outside  of  it.  A 
finite  difference  solution  to  the  Laplace  equation  is  now  desired  for 
the  point  i,j,k  taking  into  account  the  insulated  condition. 

Consider  now  the  physical  situation  which  exists  in  the  composite 
material  when  many  cells  exist  side  by  side.  The  insulated  condition 
is  a  representation  of  the  symmetry  that  exists  in  the  temperature 
field  across  the  plane.  This  symmetry  means  that  the  temperature  at 
the  point  i  +  l,j,k  which  actually  exists  as  a  point  in  the  adjacent 
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Figure  26.  Insulated  boundary  point. 

unit  cell  will  be  the  same  as  the  temperature  at  the  point  i  -  l,j,k. 
Similarly,  the  temperatures  of  the  points  i  +  2,j,k  and  i  -  2,j,k  will 
be  the  same  so  that  complete  symmetry  between  the  cells  is  achieved. 

If  a  node  is  considered  to  be  at  the  point  i  +  l,j,k  and  to  have 
the  temperature  of  point  i  -  l,j,k,  then  the  point  i,«j,k  is  seen  to 
have  six  neighbors  for  the  purpose  of  calculation.  The  general  finite 
difference  equation  is  then  applicable  here  and  can  be  written  as 
below  with  the  appropriate  equality  of  temperatures  introduced. 


(104) 


Ti,j,k=  (2Ti-l,j,k  +  Ti, j+1 ,k  +  Ti,j-l,k  +  Ti,j,k+1  + 

Ti,j,k-1)/6  (105) 

Similar  expressions  are  of  course  obtained  when  the  insulated  con¬ 
dition  occurs  on  a  plane  other  than  the  x- plane. 
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3.3-3  Interface  Points 

The  last  condition  which  must  be  satisfied  occurs  at  the  fiber- 
matrix  interfaces.  The  condition  is  that  the  temperatures  in  the 
two  regions  must  be  the  same  at  the  interface  and  that  the  heat  flux 
normal  to  the  interface  must  be  the  same  in  each  region. 

Consider  the  following  section  of  an  interface  in  which  additional 
nodal  points  have  been  added  where  the  regular  grid  lines  cross  the 
interface.  The  assumption  will  be  made  throughout  this  analysis  the 

Ax  =  Ay  • 


Figure  27.  Interface  point. 


It  must  be  kept  in  mind  that  the  z  direction  goes  into  the  paper, 
and  thus  similar  arrangements  of  nodes  lie  in  front  of  and  behind  t  e 
one  shown  in  Figure  27. 


It  is  desired  to  find  the  temperature  at  i,j,k  given  the  temp¬ 
eratures  of  all  other  nodal  points.  The  temperature  must  be  based 
solely  on  the  temperatures  of  those  points  located  in  the  same  region 
or  on  the  interface.  The  points  in  the  other  region  cannot  he  used 
because  there  is  a  change  in  conductivity  at  the  interface.  Thus, 
the  temperature  at  i,j,k  must  be  found  from  that  at  points  i  -  l,J,kJ 
i  i-1  k«  i, j ,k+l:  i,j,k-lj  aj  and  b.  This  requires  a  solution  to  the 
Laplace  equation  when  unequal  spacings  are  involved.  In  this  W> 
the  temperature  at  the  points  next  to  the  interface  can  be  found  al- 
though  no  effort  has  yet  been  made  to  satisfy  the  heat  flux  requirements 
at  the  interface.  The  temperatures  of  the  points  on  the  interface  will 
be  calculated  so  as  to  satisfy  those  requirements. 
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A  finite -difference  formulation  of  the  Laplace  equation  using 
unequal  spacings  will  now  be  found.  The  general  result  can  then  be 
used  to  find  the  temperatures  of  any  points  which  adjoin  points  on  the 
interface  and  thus  have  unequal  spacings  involved. 

Consider  the  general  arrangement  of  a  nodal  point  with  six  neighbors 
positioned  on  an  orthogonal  set  of  axes  in  Figure  28.  The  notation  is 
such  that  T  refers  to  the  temperature  at  the  specified  point  and  D 
refers  to  the  distance  between  points.  The  subscripts  L,  R,  U,  D,  F, 
and  B  refer  to  the  directions  left,  right,  up,  down,  forward,  and  backward. 


From  these,  an  expression  for  the  second  derivative  with  respect 
co  x  at  the  point  i,«j,k  can  be  obtained. 


ill!  -  l±j,  i.k  ~  dx  i-^.  ,1.1c  (108) 
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Similar  expressions  can  be  written  for  the  y  and  the  z  directions. 
These  expressions  can  now  be  inserted  into  the  Laplace  equation  to  obtain 
a  representation  of  that  equation  in  finite-difference  form. 
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This  equation  can  then  be  solved  to  obtain  an  expression  for  TQ  in 
terms  of  the  six  spacings  and  six  temperatures.  That  solution  is  as 
follows : 


To  = 
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(VDk) 
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(112) 


It  can  be  noted  that  this  equation  reduces  to  the  general  equation, 
103  when  all  of  the  spacings  are  equal.  By  use  of  this  equation,  all 
those  points  next  to  a  fiber-matrix  interface  can  now  be  calculated 
according  to  the  governing  Laplace  equation. 
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The  temperature  at  the  nodal  points  which  lie  on  the  fiber-matrix 
interface  can  be  assigned  so  as  to  satisfy  the  condition  of  equal  heat 
glux  normal  to  the  interface  in  each  region.  This  condition  can  be 
described  as  follows  where  n  is  the  direction  normal  to  the  interface: 


qI  =  qII 


(H3) 


KlJnI  ■  K 


*TlI 

ii  5^ 


(114) 


Refer  now  to  the  general  section  of  an  interface  in  Figure  29  where 
the  temperature  at  point  0  is  to  be  determined  from  those  nearest 
points  around  it.  Since  the  interface  is  in  the  shape  of  a  quarter 
circle,  the  normal  direction  is  the  same  as  the  r  direction  when  polar 
coordinates  sire  used.  Using  this  fact  and  simple  geometric  identities. 
Equation  ll4  can  be  rewritten  as  in  Equation  115- 

y  fi Tt  dx  .  ^Tt  Vvx  _  .,  /JTtt  ix  .  dTjrdyx  ,,,cv 

Ki(ir  5?  +  JT  p  =  kii(3Ti  SF  +  J,  f?>  <U5> 

It  should  be  noted  that  the  normal  direction  is  perpendicular  to 
the  z  direction  and  so  derivatives  with  respect  to  the  z  direction  do 
not  enter  into  the  problem.  The  expressions  for  the  derivatives  below 
(116  through  121)  can  now  be  inserted  into  Equation  115  to  obtain  a 
finite-difference  representation  of  that  equation.  The  derivatives 
are  defined  in  finite -difference  form  by  the  use  of  the  four  points 
right,  left,  up,  and  down  with  respect  to  point  0.  The  point  to  the 
left  is  point  3,  and  the  point  to  the  right  is  point  8.  The  point 
above  is  point  5  and  its  temperature  can  be  estimated  by  interpolation 
between  points  1  and  7*  The  temperature  of  point  6  below  point  0  can 
be  estimated  in  a  similar  manner. 
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KI[(^^)c°s(0)  +  (^)sin(e)]  = 
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11  UR 

This  finite-difference  representation  can  now  be  solved  for  TQ 
to  give  the  following  expression: 


T  T  T  T 

K^^pcos  ( 0 )  +  j=psin(e)3  +  Kjj^fpcosCe)  +  g^sin(e)] 


(123) 


By  setting  the  point  on  the  interface  to  this  value,  the  dif¬ 
ference  in  conductivities  between  the  materials  is  taken  into  account. 
Thus,  it  can  be  seen  that  the  purpose  of  putting  these  points  on  the 
interface  was  not  to  improve  the  accuracy  of  the  solution,  but  rather 
it  was  simply  to  institute  the  boundary  condition. 


Thus,  the  basic  equations  governing  every  point  in  the  unit  cell 
have  been  developed.  A  computer  can  now  be  used  to  calculate  the 
temperatures  of  all  points  in  turn.  Then  by  using  the  new  temperature 
estimates,  it  can  keep  recalculating  all  the  temperatures  until  they 
converge  to  stable  values  which  are  approximations  to  the  temperatures 
at  the  respective  points  in  the  actual  composite  material. 

3.4  SOLUTION  METHOD 

3-4.1  General  Discussion 

Two  Fortran  programs  were  written  for  the  arrangement  where 
alternate  layers  of  fibers  are  crossed  at  90-degree  angles.  They  both 
appear  in  the  appendix  to  this  report. 


The  first  case  was  where  the  heat  transfer  is  axial  to  one  layer  of 
fibers  and  transverse  to  the  adjacent  layers  of  fibers.  This  boundary 
value  problem  was  presented  earlier  in  this  section. 

The  program  models  this  problem  by  placing  a  regular  grid  on  it 
where  Ax  =  Ay  =  Az  .  The  proper  equation  is  used  for  each  point,  and 
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the  computer  succesively  calculates  the  temperature  of  each  point  while 
at  all  times  using  the  most  recent  values  of  temperatures  at  the 
neighboring  points. 

The  computer  keeps  iterating  through  the  entire  set  of  points 
until  no  temperature  of  any  point  changes  more  than  a  specified  value 
in  a  specified  number  of  iterations.  It  is  then  determined  that 
adequate  convergence  has  been  achieved.  Then  the  computer  performs  a 
numerical  integration  across  either  the  top  face  or  the  bottom  face, 
depending  on  which  is  more  accurate.  These  factors  which  affect 
accuracy  are  discussed  in  detail  in  the  following  section. 

ihe  heat  flux  is  found  according  to  the  following  equation  which 
can  be  broken  up  into  that  part  in  the  fiber  and  that  part  in  the 
matrix: 


«  ■  'Ar..<Ka»,dA 


(124) 


Q  =  Kt  Z  (77)  A  +  K  ,  £  (££)  A 
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The  value  of  AT  is  found  by  a  two  point  derivative,  or  just  the 
difference  between  the  temperatures  of  the  two  adjacent  points  at  one 
of  the  sides  of  the  unit  cell.  These  values  of  AT  are  summed  ac¬ 
cording  to  Equation  125  and  the  appropriate  values  of  Ay  and  A  A  , 
depending  on  the  grid  used,  are  inserted.  The  effective  conductivity 
is  then  found  by  Equation  126  where  AT  is  the  macro  temperature 
gradient  across  the  unit  cell  and  An  is  the  unit  cell  dimension  in 
the  direction  of  the  heat  flow. 


*eff  =  a/A(^)  '126> 

'IVo  methods  were  used  to  reduce  the  calculation  time  of  the  com¬ 
puter.  The  first  was  to  use  a  very  coarse  grid  with  large  spacings  and 
few  points  and  an  initial  guess  of  constant  temperature  to  give  a  quick 
estimate.  A  finer  grid  was  then  established  using  the  values  from  the 
coarse  grid  as  initial  guesses.  Ihe  computer  then  iterated  through  to 
the  desired  level  of  convergence.  This  downsizing  of  the  grid  was 
performed  several  times  and  resulted  in  a  much  faster  program  than  if 
the  final  grid  was  used  from  the  start  with  a  constant  temperature  initial 
guess . 

The  second  method  of  speeding  up  the  program  was  by  a  method  known 
as  over-relaxation.  Here  the  change  in  the  estimate  of  the  temperature 
at  a  point  is  calculated  according  to  those  equations  derived  previously. 
Then,  instead  of  changing  the  t.oinjvraturo  by  that  amount,  it  is  changed 
by  that  ;u:iount  plus  an  additional  fraction  of  that  amount. 


The  general  equation  (Equation  103)  would  be  altered  in  the  follow¬ 
ing  manner,  and  the  other  equations  from  Section  3*3  would  be  altered 
similarly. 


n,new 

Ai, j.k 


(Ti+l,j,k  +  Ti-l,j,k  +  Ti , j+1 ,k  +  Ti,j-l,k  + 

Ti,  j.k+l  +  Ti, 


(127) 


Ti » j  *k  A(Ti+l,j,k  +  Ti-i,j,k  +  Ti,j+l,k  +  Ti,3-l,k  + 


(128) 


Here,  A  is  known  as  the  relaxation  factor.  If  it  is  1  then 
Equation  128  reduces  to  that  of  127.  However,  if  it  is  between  1  and 
2  the  temperature,  j5k»  will  be  changed  by  a  factor  of  A  more  than 
it  would  have  otherwise,  if  a  value  of  2  or  greater  is  used  then 
numerical  oscillations  will  occur,  and  the  problem  will  not  converge. 

A  relaxation  factor  of  1.5  was  found  to  be  optimum  in  terms  of  the 
number  of  iterations  needed  to  achieve  convergence. 

3.4.2  Temperature  Field  Symmetry 

The  program  for  the  case  where  heat  flow  is  transverse  to  both 
fibers  works  in  a  similar  manner  to  the  one  described  above.  There  is, 
however,  one  additional  feature  that  can  be  used  to  reduce  the  amount 
of  computer  time  that  is  needed  to  reach  convergence.  This  feature  is 
an  additional  plane  of  symmetry  that  exists  within  the  unit  cell. 

Consider  Figure  30  which  is  the  unit  cell  with  the  midplane,  which 
is  the  plane  of  symmetry,  outlined  in  heavy  lines.  Also,  shown  are 
representative  lines  of  constant  temperature  as  they  would  appear  on 
the  front  and  left  faces  of  the  cell. 

The  temperature  has  been  arbitrarily  taken  as  +1  on  the  top 
surface  and  -1  on  the  bottom  surface.  First,  consider  point  d  on  the 
corner  of  the  cell.  It  can  be  seen  that  due  to  symmetry  the  con¬ 
ductivity  changes  above  it  exactly  counteract  the  conductivity  changes 
below  it,  and  so  it  must  be  at  a  temperature  of  zero.  The  same  reasoning 
holds  for  point  b  as  well  as  all  other  points  on  the  heavy  line  from 
d  to  b. 
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Figure  30.  Symmetry  for  transverse-transverse  heat  flow. 

Consider  now  two  points  which  lie  on  the  midplane  and  which  are 
symmetric  with  respect  to  the  line  b-d.  Points  a  and  c  are  two  such 
points.  It  can  be  seen  that  the  conductivity  effects  above  c  are  the 
same  as  those  below  a.  Similarly,  the  conductivity  effects  below  c 
arc  the  same  as  those  above  a.  Thus,  a  symmetry  exists  in  the 
temperature  field  at  these  two  points.  It  is  such  that  the  temperatures 
will  be  symmetric  with  respect  to  the  temperature  on  the  diagonal  b-d, 
and  this  is  caused  by  the  constant  temperature  lines  on  the  faces  of  the 
unit  cell.  Because  the  boundary  condition  temperatures  were  chosen  in 
such  a  way  so  as  to  make  the  temperature  on  the  diagonal  b-d  equal  to 
zero,  then  temperatures  of  symmetric  points  in  the  midplane,  such  as 
a  and  c,  will  be  negatives  of  each  other.  In  addition,  because  of 
the  symmetry  in  the  temperature  field  at  these  points  the  heat  flux 
will  be  the  same,  and  so  the  value  of  cVr/ity  will  be  the  same. 

Therefore,  only  one-half  of  the  unit  cell  needs  to  be  used  to 
obtain  a  solution  to  this  problem.  The  above  described  conditions  are 
then  required  of  the  points  on  the  midplane  so  that  a  completely  de¬ 
fined  boundary  value  problem  exists.  This  is  a  very  similar  approach 
to  that  used  to  obtain  a  reduction  of  the  unit  area  used  for  the 
staggered  array  with  parallel  fibers.  The  benefit  of  this  approach  is 
that  only  half  as  many  points  need  to  have  their  temperatures  cal¬ 
culated  for  each  interation.  Thus,  the  computer  time  needed  for  con¬ 
vergence  is  reduced  by  approximately  half. 
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There  are  a  number  of  l'.ir  tor::  which  uf  fee  t,  the  accuracy  of  the 
results  of  those  fini  t.o-di  fference  calculations.  'llu-y  include  the 
number  of  point.;-,  used,  the  plane  of  integration  for  the  heat  flux  cal¬ 
culation,  and  the  criteria  used  to  determine  when  adequate  convergence 
has  been  achieved. 

J.S.L  Number  of  Points 

ihe  number  of  ^>iut.s  used  is  a  critical  factor  affecting  the  ac¬ 
curacy  of  tile  result..  in  general,  the  result  is  more  accurate-  for  a 
larger  number  of  i-oints  because  with  more  points  the  spacing;-,  between 
points  are  smaller  and  so  the  f iuite-di fference  derivatives  more  closely 
represent  tile  true  values.  Tims,  one  would  like  to  use  a  very  large 
number  of  points.  However,  as  the  number  of  joints  increases,  the 
computer  time  needed  increases  sharply,  ;uid  eventually  the  result  could 
become  worse  due  to  roundoff  errors.  All  of  these  differences  must  be 
added  in  the  integration  process  with  roundoff  occurring  after  each 
operation,  and  with  tin-  number  of  points  being  increased,  the  number  of 
calculations  is  Increased. 

A  finite-difference  program  ■..•as.  written  for  the  roetiuigular  array 
of  feet. ion  l  i  to  study  the  effects  of  the  various,  pjiramotors  on  the 
accuracy  of  the  result..  A  very  good  result  had  been  obtained  by  other 
methods,  to  which  the  f  ini  te-di fference  resul  ts,  could  tie  compared.  As 
;m  example,  all  paivimetors  were  kept,  (lie  s;uue  except  the  munber  of 
di  vis,  ions  along  each  side  of  the  cell  was.  varied. 

TAhlJ-:  0.  K.KI4VTU  OK  NUMHKK  OK  htVIOlONf- 

Number  of  Pi  visions  ’VffAmatrix  lnterations 
('er  f,  ide 


10 

l.i-‘p!lo 

j:x> 

ps 

l.o^lo 

1000 

so 

1.0U50 

j.-'oo 

100 

i.tvW 

7800 

Accurate  Value  Using  Other  Method *Vff  Amatri x  -  l»b5^0 

'Ihese  data  demonstrate  the  decrease  In  accuracy  for  a  e-mail  munber 
of  divisions  and  for  a  large  number  ot’  divisions.  The  value  for  a  high 
number  of  divisions  can,  however,  be  greatly  improved  by  requiring  a 
better  tolerance  for  eacli  of  the  individual  temperatures .  Hy  iterating 
through  more  times  to  obtain  more  significant  figures  in  each  of  the 
temperature  values,  the  effect  of  roundoff  errors  can  be  reduced  sub¬ 
stantially.  liiis  will  of  course  increase  the  computing  time  even  more 
than  has  already  accumulated  because  of  the  large  munber  of  points.  For 
these  reasons  and  because  some  of  the  figures  above  are  correct  to 
four  significant  figures  already,  it  does  not.  seem  necessary  to  use  an 
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extremely  large  number  of  points.  For  the  3-dimensional  problem  considered 
here,  28  points  per  side  were  used  along  the  base  and  the  number  of  points 
used  along  the  vertical  dimension  was  determined  by  the  dimensions  of  the 
specific  problem. 

3.5.2  Convergence  Criterion 


The  criterion  that  is  used  to  determine  when  the  computer  should  stop 
iterating  the  temperatures  of  the  points  has  a  significant  effect  on  the 
accuracy  of  the  result.  The  program  was  designed  to  stop  when  every 
temperature  changed  no  more  than  a  specified  amount  in  20  iterations.  A 
figure  of  20  was  chosen  to  guarantee  that  convergence  really  had  been 
reached  and  that  all  the  temperatures  had  stopped  fluctuating  more  than 
the  specified  amount.  That  specified  amount  will  now  be  referred  to  as 
the  tolerance . 


If  the  tolerance  is  too  large,  then  when  the  program  stops  there 
will  not  be  enough  significant  figures  in  the  temperature  values  to 
survive  the  effects  of  roundoff  in  the  numerical  integration.  If  the 
tolerance  is  too  small,  then  the  time  that  the  program  runs  may  become 
excessive,  'lhe  following  table  shows  the  effect  of  tolerance  on  the 
accuracy  of  the  result  and  the  number  of  iterations  required  to  achieve 
it.  'Ihis  data  is  from  the  2-dimensional  finite-difference  program  which 
models  the  rectangular  array  of  Section  ! ! . 

TAB  IF.  4.  EFFECTS  01  TOiFKANCE 


Tolerance 

^eff/^matrix 

Iterations 

.00! 

1.6489 

580 

.0001 

1.051b 

1000 

.00005 

1.0523 

1160 

.00001 

1.6529 

1420 

Accurate  Value  Using  Other  Method - K^ff/^matrix  =  1*6520 

For  the  results  presented  in  the  appendix  to  this  report,  a  value 
of  0.0001  was  used.  In  the  table  above,  this  provided  a  result  accurate 
to  four  significant  figures.  Further,  another  order  of  magnitude 
improvement  in  the  tolerance  did  not  provide  any  more  significant 
figures.  Since  three  significant  figures  is  adequate,  a  value  of  0.0001 
was  more  than  satisfactory,  and  approximately  three-figure  accuracy  was 
obtained  for  all  the  results  presented  in  the  appendix. 

3.5*3  I ntegration  Plane 

'lhe  other  factor  that,  can  affect  the  accuracy  of  the  result  is 
the  choice  of  which  plane  to  use  for  the  numerical  integration  to  find 
the  heat  flux,  'lhe  surface  planes  of  the  unit  cell  are  the  best  planes 
to  use  because  either  she  fiber-matrix  interface  is  not  present  there  or 
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it  intersects  the  plane  at  a  right  angle.  Either  way,  these  planes  provide 
for  the  simplest  numerical  integration. 

The  accuracy  of  the  result  deteriorates  as  the  volume  ratio  in¬ 
creases  and  as  the  ratio  of  fiber  conductivity  to  matrix  conductivity 
is  increased.  Consider  the  following  2-dimensional  finite-difference 
representation  which  will  show  the  same  effects  that  occur  with  a  3- 
dimensional  problem; 


T  =  0 

Figure  31>  high  volume  ratio  array. 

In  this  case  a  high  volume  ratio  is  used,  and  only  two  nodal 
points  exist  on  a  line  from  the  fiber  to  the  edge  of  the  unit  cell. 

Assume  also  that  a  high  conductivity  ratio  lias  been  used.  Because 
of  these  conditions,  most  of  the  heat  flow  along  the  top  surface  will 
occur  in  those  nodal  points  on  the  left.  On  the  bottom  surface  though } 
the  heat  flow  will  be  spread  throughout  the  fiber,  and  thus  it  will 
occur  through  a  much  highei  number  of  points.  For  this  case,  integration 
across  the  bottom  surface  will  give  a  more  accurate  answer  than  across 
the  top  surface.  The  result  from  an  integration  across  the  top  will  be 
less  accurate  because  there  are  not  as  many  points  to  give  a  description 
of  the  temperature  field  in  that  region  that  will  contribute  the  most 
to  the  heat  flux. 

Consider  now  the  case  where  there  are  more  nodes  between  the  fiber 
and  the  top  edge  in  Figure  $2.  Again,  assuming  a  high  conductivity 
ratio,  the  heat  flow  will  now  be  spread  out  through  the  top  surface 
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much  more  uniformly  than  was  the  case  before.  Further,  the  integration 
across  the  bottom  will  be  worse.  This  is  because  most  of  the  resistance 
to  the  heat  flow  is  occuring  in  the  matrix,  so  most  of  the  temperature 
drop  will  occur  at  those  points.  Thus,  the  fiber  will  be  very  close  to 
a  constant  temperature.  When  the  differences  are  found  between  the 
temperatures  on  the  bottom  plane  and  those  on  the  adjacent  plane  above 
it,  the  differences  will  be  very  small  and  so  roundoff  errors  will 
reduce  the  accuracy.  Therefore,  for  this  case  the  integration  is  best 
performed  along  the  top  edge. 


Figure  32.  Low  volume  ratio  array. 

It  should  be  noted  that  roundoff  errors  are  less  when  the  case  of 
only  a  few  points  between  the  fiber  and  edge  is  considered.  This  is 
because  a  higher  portion  of  the  resistance  to  heat  flow  is  then  located 
in  the  fiber  region.  Thus,  there  will  be  more  of  a  temperature  drop 
in  the  fiber  than  when  almost  all  of  the  resistance  is  located  in  the 
matrix.  Thus,  the  top  surface  should  be  used  when  there  are  a  number 
of  nodes  between  the  fiber  and  edge.  Otherwise  the  bottom  surface 
should  be  used. 

Consider  now  the  case  where  heat  flow  is  transverse  to  both  fibers. 
The  unit  cell  that  was  modelled  by  the  computer  is  shown  below. 
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Figure  33-  0°-90°  Unit  cell  for  heat  flow 

transverse  to  both  fibers. 


This  problem  is  much  like  the  2-dimensional  case  described  above. 
Thus,  integration  should  be  performed  according  to  the  considerations 
discussed  for  it.  It  should  be  noted  though  that  the  two  integration 
schemes  will  give  better  agreement  as  the  conductivity  ratio  is 
reduced. 


Consider  now  the  case  of  heat  flow  axial  to  one  fiber  and  transverse 
to  the  other.  The  axial  fiber  will  carry  most  of  the  heat  flow  for 
high  conductivity  ratios.  Because  of  its  orientation,  the  temperature 
drop  from  front  to  back  will  be  quite  uniform  along  its  length.  Also 
because  of  this  fiber's  dominant  effect  on  the  heat  transfer,  the 
temperature  drop  through  the  entire  cell  will  tend  to  be  more  uniform 
than  was  the  case  for  the  previous  problem.  Thus,  the  roundoff  errors 
on  the  front  face  will  not  be  as  severe  as  they  were  across  the  bottom 
face  of  the  previous  problem  since  the  temperature  differences  between 
adjacent  points  are  much  greater.  In  addition,  with  this  arrangement 
most  of  the  heat  flux  is  occurring  across  the  axial  fiber  which  is  a 
much  larger  region  than  was  the  case  for  the  previous  problem,  'fliis 
would  tend  to  indicate  that  integration  across  the  back  plane  would 
provide  an  accurate  result. 

The  actual  results  showed  that  both  methods  were  ac evirate  and  that 
they  agreed  to  three  significant  figures  for  all  the  results  presented 
in  the  appendix  to  this  report.  This  showed  that  the  effects  of  roundoff 
errors  on  the  front  face  were  just  as  small  as  the  effects  of  heat 
transfer  through  a  small  number  of  points  on  the  back  face. 
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3-6  RESULTS 


The  results  for  the  0°-90°  array  are  presented  in  graphical  form  in 
Appendix  A.  These  graphs  present  the  effective  conductivities  for  aspect 
ratios  of  1  and  2  (layer  spacing  to  fiber  spacing).  Further,  the  con¬ 
ductivity  ratio  is  varied  from  0.1  to  100,  and  the  volume  ratio  takes 
on  the  values  0.4,  0.5,  0.6,  O.65  and  0.70.  The  spacing  between  fibers 
in  adjacent  layers  is  taken  to  be  the  same  throughout  the  material. 

Of  the  three  principal  conductivities  in  these  materials,  two  will 
always  be  the  same.  These  are  the  two  in  which  heat  flow  is  transverse 
to  one  layer  of  fibers  and  parallel  to  the  adjacent  layers  of  fibers. 
Parametric  curves  for  these  conditions  are  presented  in  Figures  A-59 
and  A-60.  The  third  principal  conductivity  is  such  that  the  heat  flow 
is  transverse  to  all  fibers,  and  these  results  are  presented  in 
Figures  A-57  and  A-58.  The  results  from  which  all  of  these  curves  were 
drawn  were  accurate  to  at  least  three  significant  figures. 

For  the  cases  where  heat  flow  is  transverse  to  both  layers  of 
fibers,  it  can  be  seen  that  the  curves  level  off  at  high  conductivity 
ratios.  This  is  because  all  of  the  heat  must  pass  through  a  portion 
of  the  matrix,  and  at  high  conductivity  ratios  almost  all  of  the 
resistance  to  heat  flow  is  in  the  matrix.  Thus,  decreasing  the  resist¬ 
ance  of  the  fiber  still  further  affects  the  total  resistance  that  the 
heat  fiasses  through  very  little,  and  the  effective  conductivity  is 
increased  very  li ttle . 

Consider  now  the  case  where  heat  flow  is  axial  to  one  layer  of 
fibers  and  transverse  to  the  other.  There  the  axial  fibers  provide 
a  very  low  resistance  piath  for  the  heat  flow  when  the  conductivity  ratio 
is  high.  Most  of  the  heat  flow  will  be  within  the  axial  fiber,  and  as 
the  conductivity  of  that  fiber  is  increased,  the  overall  effective 
conductivity  will  Increase  by  the  same  order  of  magnitude.  Thus,  the 
curves  for  the  axial- transverse  case  are  steeply  sloped  at  high  con¬ 
ductivity  ratios  and  will  approach  linear  lines  as  the  additional  con¬ 
tribution  of  the  transverse  fiber  becomes  negligible. 

Consider  the  graph  of  transverse- transverse  heat  flow  for  an  aspect 
ratio  of  1.0.  If  alternate  layers  of  fibers  wore  turned  90  degrees  the 
45  degree  case  for  the  parallel  fibers  in  a  rectangular  array  would  be 
obtained.  [n  this  way  the  0°-90°  array  is  similar  to  the  rectangular 
array  of  45  degrees . 

By  comparing  one  curves  of  these  two  cases  it  can  be  seen  that  they 
are  identical  for  volume  ratios  less  than  or  equal  to  0.40.  This  shows 
that  at  low  volume  ratios  the  orientation  of  the  fibers  does  not  affect 
the  effective  conductivity  significantly .  As  the  volume  ratio  is  in¬ 
creased,  the  curves  for  the  0°-90°  case  fall  further  and  further  below 
the  corresponding  curves  for  the  rectangular  array.  'Jhis  is  because  the 
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heat  conduction  paths  in  the  0°-90°  array  are  not  as  straight,  and  the 
heat  is  forced  to  pass  through  a  greater  amount  of  resistive  matrix 
material.  When  the  fibers  are  parallel,  it  can  pass  much  more  directly 
from  one  fiber  to  the  next  and  pass  through  a  minimum  amount  of  matrix 
material. 

For  the  0°-90°  array  with  an  aspect  ratio  of  2,  alternate  layers 
can  be  turned  to  produce  a  rectangular  array  with  a  =  63.4°.  Urns, 
these  two  materials  are  much  alike.  As  is  expected,  the  curves  for  these 
materials  are  identical  for  volume  ratios  less  than  0.39  which  is  the 
highest  value  that  is  physically  possible.  Thus,  for  low  volume  ratios, 
the  data  for  parallel  fibers  can  be  used  for  other  orientations  as  the 
effective  conductivity  is  almost  entirely  independent  of  the  orientation. 


Figure  A- 2  Parallel  Fiber  (Rectangular  Array)  Volume  Ratio 


Figure  A- 3  Parallel  Fiber  'Rectangular  Array)  Volume  Ratio 


Figure  A-4  Parallel  Fiber  (Rectangular  Array)  Volume  Ratio 


Figure  A- 9  Parallel  Fiber  (Rectangular  Array)  Volume  Ratio 
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"  fiber  '  "  matrix 

Figure  A- 12  Parallel  Fiber  (Rectangular  Array)  Volume  Ratio 


Figure  A- 14  Parallel  Fiber  (Rectangular  Array)  Volume  Ratio 


Figure  A-l8  Parallel  Fiber  Rectangular  Array)  Volume  Ratio 
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Figure  A-21  Parallel  Tiber  (Staggered  Array)  Volume  Ratio 
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Figure  A-24  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 


Ratio 


i-30  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 
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Figure  A- 33  Parallel  Fiber  (staggered  Array)  Volume  Ratio 


Figure 


fiber  '  ^  motrix 

Figure  A- 38  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 


Figure  A-39  Parallel  Fiber  (staggered  Array)  Volume  Ratio  =  .60 


Figure  A-40  Parallel  .fiber  ( Staggered  Array)  Volume  Ratio 
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Figure  A-4l  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 


l-42  Parallel  Fiber  (Staggered  Array)  Volume 
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Parallel  Fiber  (Staggered  Array)  Volume  Ratio 
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Figure  A-47  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 


Figure  A-48  Parallel  Fiber  (Staggered  Array)  Volume  Patio 


Figure  A-50  Parallel  Fiber  ( Staggered  Array)  Volume  Ratio 


Figure  A-51  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 


Figure  A-52  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 


Figure  A- 53  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 


{-56  Parallel  Fiber  (Staggered  Array)  Volume  Ratio 


Figure  A-57  0°  -  90°  Array  (Transverse  -  Transverse  Array)  Aspect  Ratio 


Figure  A- 58  0°  -  90°  Array  'Transverse  -  Transverse  Array)  Aspect  Ratio 


Figure  A-6l  0°  -  90°  Array  -  Volume  Ratio  =  .200,  Aspect  Ratio 


A-62  0°  -  90°  Array  -  Volume  Ratio  =  .250,  Aspect  Ratio 


Volume  Ratio  =  .300,  Aspect  Ratio 


A- 65  0°  -  90°  Array 


Array  -  Volume  Ratio  =  .500,  Aspect  Ratio 


Axial  -  Transverse 
Transverse  -  Transverse 
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Figure  A-68  0°  -  90°  Array  -  Volume  Ratio  =  .600,  Aspect  Ratio 


Figure  A-69 


Figure  A- 70  0°  -  90°  Array  -  Volume  Ratio  =  .700,  Aspect  Ratio 
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APPENDIX  B 


PARALLEL  FIBER 
PROGRAM  METHOD 


The  computer  programs  which  were  written  to  calculate  the  effective 
conductivity  for  both  the  staggered  and  retangular  arrays  follow  the 
methods  of  Sections  2.2  through  2.4.  In  the  manner  discussed,  a  series 
of  M  equations  with  M  unknowns  are  produced  where  M  is  the  number  of 
discrete  points  which  are  used  to  satisfy  the  boundary  conditions. 

These  equations  are  solved  for  the  unknowns  which  are  then  put  into  the 
appropriate  equations  to  give  the  effective  conductivity  for  the  material 
under  study. 

The  coefficients  are  found  by  the  use  of  the  appropriate  boundary 
condition  equations.  The  M  equations  can  then  be  written  using  matrix 
notation  in  the  following  manner: 

^MxM^Mxl  =  (C^Mxl 

Here,  A  is  the  M  by  M  coefficient  matrix,  E  is  the  solution  matrix 
(actually  it  is  just  a  vector),  and  C  is  the  matrix  (also  a  vector) 
containing  the  constants  from  the  boundary  conditions. 

The  program  employs  a  gauss-elimination  technique  whereby  rows  are 
multiplied  through  by  constants  and  subtracted  in  such  a  way  so  as  to 
produce  an  upper  triangular  matrix.  The  unknowns  can  then  b<  found  one 
at  a  time  by  starting  at  the  bottom  and  back-substituting  the  unknowns 
into  the  equations  as  they  are  found.  Reference  7  gives  a  very  good 
description  of  this  method. 

A  standard  gauss -elimination  procedure  is  not  extremely  accurate 
because  of  roundoff  errors  within  the  computer.  This  can  be  greatly 
improved  by  a  procedure  called  pivoting  where  rows  are  interchanged 
during  the  elimination  process  to  put  the  largest  terms  in  each  column 
on  the  main  diagonal.  The  programs  used  here  went  one  step  further 
and  included  a  double  pivoting  procedure  where  throughout  the  elimination 
the  largest  terms  in  the  entire  array  are  put  in  the  key  positions  on 
the  diagonal.  This  means  that  columns  are  interchanged  as  well  as  rows. 
This  procedure  requires  an  intricate  bookkeeping  scheme  within  the 
program,  but  it  decreases  the  effects  of  roundoff  erroic  still  further 
beyond  what  the  single  pivoting  can  do.  The  details  of  these  methods 
will  not  be  discussed  here  as  they  are  covered  in  many  references 
(again,  Reference  7  is  quite  good). 

Another  method  used  to  decrease  the  effects  of  roundoff  errors  in 
these  programs  was  the  use  of  double  precision  logic.  This  doubled 
the  number  of  digits  which  were  carried  throughout  the  calculations 
within  the  computer  from  8  to  16  and  forced  the  err^r  to  occur  in  the 
last  of  these  digits.  \ 

In  the  manner  described  above,  every  possible  technique  has  been 
used  to  limit  the  error  due  to  roundoff  after  each  calculation  by  the 
computer.  However,  there  is  another  type  of  error  which  is  due  to  the 
fact  that  the  power  series  which  describe  the  temperature  field,  are 
chopped  off  after  a  finite  number  of  terms. 


Theoretically,  it  is  expected  that  this  error  in  the  solution 
would  be  reduced  when  more  and  more  terms  in  the  power  series  are 
taken  into  consideration.  This  would  be  true  if  an  infinite  number  of 
digits  were  carried  through  all  the  calculations.  However,  given  a 
fixed  number  of  digits,  a  point  is  reached  after  which  the  increase  in 
accuracy  due  to  added  terms  in  the  series  is  more  than  offset  by  in¬ 
creased  error  due  to  roundoff.  When  several  terms  are  added  past  this 
point,  the  roundoff  error  increases  rapidly  to  the  point  where  the  method 
becomes  unstable  and  the  result  becomes  meaningless. 

The  error  in  the  solution  is  high  when  a  small  number  of  discrete 
matching  points  are  used  due  to  the  chop-off  of  the  power  series.  It 
is  again  high  when  a  large  number  of  points  are  used  due  to  roundoff, 
and  so  there  is  an  optimum  number  of  points  that  should  be  used.  This 
optimum  number  was  found  to  be  15  points  for  the  rectangular  array  and 
approximately  that  number  for  the  staggered  array.  The  locations  of 
these  points  are  discussed  in  Section  II.  The  number  of  points  for  the 
staggered  array  varies  somewhat  due  to  the  intricacies  of  the  matching 
line. 


By  using  this  method  of  exactly  satisfying  the  boundary  conditions 
at  discrete  points,  it  was  hoped  that  they  are  closely  satisfied  at  the 
points  between  the  discrete  points.  As  an  example  of  the  accuracy 
obtained,  consider  the  rectangular  array.  Along  the  top  edge  the 
temperature  value  deviated  from  1  by  only  ±  0.0001,  at  most,  and  the 
resulting  heat  flux  distribution  along  the  right  side  was  at  most  ±  0.01 

percent  of  the  average  heat  flux  along  the  top  and  bottom  edges.  For 

the  staggered  array  the  deviation  from  exact  symmetry  was  on  the  order 
of  ±  0.01  percent  of  the  values  for  both  the  temperature  and  dT/dy 
symmetry. 

At  higher  values  of  (3  (fiber  conductivity/matric  conductivity)  the 
error  in  the  effective  conductivity  becomes  greater.  In  addition,  the 
error  is  increased  at  higher  values  of  the  volume  ratio.  At  the  highest 

possible  volume  ratios  and  a  p  of  100,  the  value  of  K^pp/Kma-tm'  x  is 

accurate  to  three  significant  figures.  Further,  this  increases  to 
more  than  six  significant  figures  for  the  lower  values  of  p  and  volume 
ratio.  This  conclusion  is  arrived  at  by  comparing  the  data  with  that 
from  a  routine  which  uses  finite  difference  methods  (discussed  in 
Section  III)  and  by  studying  the  effects  of  changes  in  the  number  of 
matching  points.  This  conclusion  was  further  strengthened  by  comparing 
certain  data  between  the  staggered  and  rectangular  arrays  which  happened 
to  be  coincident.  Even  though  the  two  methods  of  solution  are  distinctly 
different,  the  very  good  agreement  between  results  for  certain  geometries 
added  confidence  to  the  results  of  both  methods.  These  various 
relationships  are  discussed  in  Section  II. 
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